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We  discuss  (2+l)-dimensional  gapless  surface  theories  of  bulk  (3+l)-dimensional  topological 
phases,  such  as  the  BF  theory  at  level  K,  and  its  generalization.  In  particular,  we  put  these  theories 
on  a  flat  (2+1)  dimensional  torus  T3  parameterized  by  its  modular  parameters,  and  compute  the 
partition  functions  obeying  various  twisted  boundary  conditions.  We  show  the  partition  functions 
are  transformed  into  each  other  under  SL( 3,  Z)  modular  transformations,  and  furthermore  establish 
the  bulk-boundary  correspondence  in  (3+1)  dimensions  by  matching  the  modular  S  and  T  matrices 
computed  from  the  boundary  field  theories  with  those  computed  in  the  bulk.  We  also  propose  the 
three-loop  braiding  statistics  can  be  studied  by  constructing  the  modular  S  and  T  matrices  from 
an  appropriate  boundary  field  theory. 
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The  bulk-boundary  correspondence  is  one  of  the  most 
salient  features  of  topologically  ordered  phases  of  mat¬ 
ter.  In  topologically  ordered  states  in  (2+1)  dimensions 
[(2+l)d],  all  essential  topological  properties  in  their  bulk 
can  be  derived  and  understood  from  their  edge  theo¬ 
ries,  such  as  quantized  transport  properties,  properties  of 
bulk  quasiparticles  (fractional  charge  and  braiding  statis¬ 
tics  thereof),  and  the  topological  entanglement  entropy, 
etc. 1-8  Edge  or  surface  theories  also  play  an  important 
role  in  symmetry-protected  and  symmetry-enriched  topo¬ 
logical  phases. 9-15 

The  purpose  of  this  paper  is  to  study  the  bulk¬ 
boundary  correspondence  in  the  simplest  (3+l)d  topo¬ 
logical  field  theory,  the  BF  topological  field  theory16-22, 
and  its  generalizations.  The  BF  theory  describes,  among 
others,  the  long  wave-length  limit  of  BCS  superconduc¬ 
tors,  and  the  deconfined  phase  of  the  Zk  gauge  theory. 
It  is  also  relevant  to  the  hydrodynamic  description  of 
(3+l)d  symmetry-protected  topological  (SPT)  phases  in¬ 
cluding  topological  insulators  and  related  systems.23-31 

To  put  our  purpose  in  the  proper  context,  let  us  give 
a  brief  overview  of  the  bulk-boundary  correspondence  in 
(2+l)d  topologically  ordered  phases.  For  (2+l)d  topo¬ 
logical  phases,  bulk  topological  phases  can  be  character¬ 
ized  by  the  modular  S  and  T  matrices.  The  S  and  T 
transformations  generate  the  basis  transformation  in  the 
space  of  degenerate  ground  states,  which  appear  when 
the  system  is  put  on  a  spatial  two-dimensional  torus. 
Combined  together,  the  S  and  T  transformations  form 
the  group  SL{ 2,Z),  the  mapping  class  group  of  the  two- 
dimensional  torus  T2.  Their  geometric  meanings  are  the 
90°  rotation  and  Delrn  twist  defined  on  the  torus,  respec¬ 
tively.  In  the  basis  in  which  the  T  matrix  is  diagonal  (the 
so-called  quasi-particle  basis) ,  the  diagonal  entries  of  the 
T  matrix  encode  the  information  on  the  topological  spin 
of  quasi-particles.  On  the  other  hand,  the  S  matrix  con¬ 
tains  the  information  of  the  braiding  and  fusion.  For  an 
Abelian  topological  phase,  the  elements  of  the  S  matrix 
are  given  by  braiding  phases  between  quasiparticles,  up 
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to  an  over  all  normalization  factor  1/2?,  where  2?  is  the 
total  quantum  dimension. 

On  the  other  hand,  at  their  boundary  (edge),  gapless 
boundary  excitations  supported  by  a  (2+l)d  topologi¬ 
cal  phase  can  be  described  by  a  (l+l)d  conformal  field 
theory  (CFT).32  There  is  one-to-one  correspondence  be¬ 
tween  quasi-particle  excitations  in  the  bulk  and  primary 
fields  living  on  the  edge.  On  the  (l+l)d  spacetime  torus, 
one  can  form  the  character  \j  (r)  from  the  tower  of  states 
built  upon  a  primary  field  Oy. 

Xi(r)  =  Trw,  [e2-^-2-^]  (1) 

where  Hq  and  Pq  are  the  Hamiltonian  and  the  mo¬ 
mentum  operators,  respectively,  the  complex  parameter 
t  =  Ti  +  1T2  is  the  modular  parameter  parameterizing 
the  spacetime  torus,  and  the  trace  is  taken  over  all  states 
in  the  Hilbert  space  Hj  that  is  built  upon  the  highest 
weight  state  associated  with  the  primary  field  Oj.  The 
characters  Xj  transform  into  each  other  under  the  mod¬ 
ular  transformations  of  the  spacetime  torus.  Under  the 
modular  T  and  S  transformations,  the  characters  Xj(T) 
transform  as 

Xj  (r  +  1)  =  e2^-^Xj(r)  =  TnXj(r ), 

Xj(~l /r)  =  '%2sjj'Xj'(T),  (2) 
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where  the  matrices  T  and  S  represent  the  action  of  the 
T  and  S  modular  transformations  on  the  characters,  re¬ 
spectively.  The  matrix  T  is  a  diagonal  matrix  and  in¬ 
cludes  the  conformal  dimension  h  for  each  character  and 
the  central  charge  c  for  the  CFT.  The  T  and  S  matri¬ 
ces  for  the  characters  in  the  edge  theory  have  the  direct 
correspondence  (and  are  essentially  identical)  to  the  the 
T  and  S  matrices  defined  for  the  corresponding  (2  +  l)d 
bulk  topological  theory. 

Coming  back  to  our  main  focus,  i.e. ,  (3  +  l)d  topolog¬ 
ically  ordered  phases,  the  bulk  topological  system  can 
be  defined  on  a  spatial  torus  T3  (while  other  choices 
are  of  course  possible).  The  mapping  class  group  of  the 
three-dimensional  torus  is  SL{ 3,  Z),  and,  as  in  the  case  of 
(2+l)d,  is  also  generated  by  two  transformations,  which 
we  also  call  the  modular  S  and  T  transformations  (see 
Sec.  II A  for  details).  For  (3  +  l)d  topological  phases 
defined  on  a  spatial  torus  T3,  S  and  T  matrices  can  be 
introduced  to  describe  the  basis  transformation  of  degen¬ 
erate  ground  states.  As  in  (2+l)d,  the  S  and  T  matrices 
encode  the  topological  data  of  the  bulk  topological  phase, 
such  as  the  braiding  and  spin  statistics  of  excitations.33,34 
In  (3+l)d,  the  exchange  statistics  of  particles  has  to  be 
either  fermionic  or  bosonic.  On  the  other  hand,  a  particle 
and  a  loop-like  excitation,  or  two  loop-like  excitations  in 
the  presence  of  an  additional  background  loop,  can  have 
non-trivial  braiding  and  can  obey  non-trivial  statistics. 
For  the  Abelian  topological  phase  described  by  the  BF 
theory,  the  S  matrix  describes  the  braiding  phase  be¬ 
tween  particle  and  loop  excitations,  while  the  T  matrix 


has  the  physical  meaning  of  a  (3  +  l)d  analogue  of  topo¬ 
logical  spins.33  It  has  been  also  proposed  that  there  exist 
(3+l)d  topological  phases  that  are  characterized  by  their 
three-loop  braiding  statistics. 34-40 

We  will  demonstrate  that  these  results,  obtained  and 
discussed  previously  from  the  bulk  point  of  view,  can 
be  obtained  solely  from  gapless  boundary  field  theo¬ 
ries.  More  specifically,  taking  various  examples  of  (3+l)d 
topologically  ordered  phases  and  their  surface  states, 
which  we  put  on  the  (2  +  l)d  spacetime  torus  T3,  we 
compute  the  modular  S  and  T  matrices  explicitly,  and 
show  that  they  agree  with  the  S  and  T  matrices  obtained 
from  the  bulk  considerations.  We  thereby  establish  the 
bulk-boundary  correspondence  in  these  (3+l)d  topolog¬ 
ically  ordered  phases.  Along  the  course,  we  also  propose 
a  bulk  continuum  field  theory  which  realizes  non-trivial 
three-loop  braiding  statistics. 

N.B.  Our  strategy  adopted  in  this  paper  is  to  utilize 
boundary  field  theories  to  learn  about  bulk  excitations 
in  (3+l)d  topological  phases,  by  establishing  a  bulk¬ 
boundary  correspondence.  One  should  however  bear  in 
mind  that  boundaries  may  have  more  “life”  than  their 
corresponding  bulk,  in  that  a  given  bulk  topological 
phase  can  be  consistently  terminated  by  more  than  one 
boundary  theory.  Therefore,  it  would  be  more  appropri¬ 
ate  to  consider  a  “stable  equivalent  class”  of  boundary 
theories  for  a  given  bulk  theory.  (See,  for  example,  Ref. 
41.)  Nevertheless,  one  can  expect  universal  topological 
properties  of  the  bulk  theories  may  be  extracted  from  any 
boundary  theory  which  consistently  terminates  the  bulk. 


A.  Outline  of  the  paper 

The  rest  of  the  paper  is  organized  as  follows. 

In  Sec.  II,  we  consider  the  compactified  free  boson  the¬ 
ory  in  (2+l)d  defined  on  the  3d  flat  torus  T3  is  computed. 
This  (2+l)d  theory  is  not  necessarily  tied  to  a  particular 
(3+l)d  bulk  topological  order  but  serves  as  a  warm  up 
for  later  sections.  We  will  show  its  partition  function  is 
invariant  under  SL(3,  Z). 

In  Sec.  Ill,  the  surface  theory  of  the  (3+l)d  BF  the¬ 
ory  at  level  K  is  studied.  This  theory  can  be  subjected  to 
twisted  boundary  conditions,  which  are  induced  by  intro¬ 
ducing  quasi-particles  in  the  (3+l)d  bulk.  We  will  show 
that  the  partition  functions  with  different  boundary  con¬ 
ditions  are  transformed  into  each  other  under  SL( 3,Z), 
and  form  a  representation  of  SL( 3,Z).  The  extracted 
S  and  T  matrices  agree  with  the  known  result.33  We 
will  also  compute  the  thermal  entropy  in  Sec.  HID,  and 
show  that  there  is  a  constant  negative  contribution  to 
the  entropy.  This  contribution  to  the  boundary  thermal 
entropy  is  expected  to  capture  the  topological  entangle¬ 
ment  entropy  defined  in  the  corresponding  (3+l)d  bulk. 

In  Sec.  IV,  we  introduce  an  additional  term,  the  ax- 
ion  term  or  the  theta  term,  to  the  (3+l)d  BF  theory. 
The  theta  angle  has  a  texture  (spatial  inhomogeneity) 
and  affects  the  boundary  theory  by  twisting  the  quantum 
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numbers.  Being  static,  the  texture  in  the  theta  angle  is 
interpreted  as  a  topological  defect,  and  we  will  show  that 
the  introduction  of  the  defect  makes  the  surface  theory 
non- modular  invariant,  in  the  sense  that  the  action  of 
modular  transformations  is  not  closed  within  the  space 
of  the  partition  functions. 

This  BF  theory  with  the  theta  term  motivates  us  to 
consider  yet  another  theory  in  Sec.  V,  which  can  be  con¬ 
structed  by  coupling  two  copies  of  the  BF  theory.  Com¬ 
pared  with  the  the  defect  system  (the  BF  theory  with 
the  theta  term)  discussed  in  Sec.  IV,  in  the  coupled  sys¬ 
tem,  each  copy  can  be  interpreted  as  playing  a  role  of  a 
defect  to  the  other.  In  this  system,  however,  there  is  no 
externally  imposed  texture.  We  propose  this  continuum 
bulk  held  theory  realizes  three-loop  braiding  statistics 
discussed  previously.35,37-39  On  the  surface,  we  consider 
two  copies  of  the  BF  surface  theories,  which  are  coupled 
together  in  their  zero  mode  sectors.  We  will  show  that, 
by  computing  the  modular  S  and  T  matrices  explicitly, 
this  system  exhibits  three-loop  braiding  statistics. 

Finally,  we  conclude  in  Sec.  VI. 


II.  THE  COMPACTIFIED  FREE  BOSON  IN 

(2+  1)D 

The  compactihed  real  scalar  theory  in  (2  +  l)d  is  de¬ 
scribed  by  the  Lagrangian  density 

C  =  Wy  [(^ } 2  -  ^  (O,0)2]  >  (3) 

where,  for  now,  the  spacetime  is  the  “canonical”  hat  torus 
T3  parameterized  by  (t,  x,  y).  (We  will  consider,  momen¬ 
tarily  in  Sec.  II B,  a  generic  torus  parameterized  by  six 
modular  parameters.)  The  boson  held  obeys  the  com- 
pactihcation  condition  on  a  circle  of  radius  r,  i.e., 

(/>=</)  +  27T1'.  (4) 

This  model  can  be  exactly  solved  and  is  dual  to  the  com¬ 
pact  XJ{  1)  gauge  theory.  Under  the  duality,  the  boson 
held  (f>  is  related  to  the  [7(1)  gauge  held  a M  by 

e^xdvax  o  9^,  o  d^d^.  (5) 

Furthermore,  quantized  vortices  on  the  boson  side  are 
dual  to  quantized  charges  in  the  [7(1)  gauge  theory.  For 
the  compact  [7(1)  gauge  held  theory,  the  monopole  (in- 
stanton)  proliferation  leads  to  a  conhning  phase  and  this 
process  on  the  scalar  boson  side  corresponds  to  adding  a 
cos {(f>)  term.42,43  This  process  breaks  the  [7(1)  symmetry 
in  the  compact  boson  theory,  and  the  particle  number  is 
not  conserved  anymore.  If  we  prohibit  the  monopoles,  on 
the  other  hand,  the  Abelian  [7(1)  gauge  theory  is  stably 
gapless. 

The  free  boson  theory  can  be  canonically  quantized: 
The  corresponding  Hamiltonian  is 

2-kR!  2t ri?2 

H  =  J  dx  J  dy  [n2  +  ( dx<f> )2  +  [dy^f]  ,  (6) 
0  0 


where  x  and  y  are  periodic  with  radius  2nRi  and  27ri?2, 
respectively,  and  the  canonical  momentum  is  (r  :=  ( x1  y)) 

H(r)  :=  dt<f>( r).  (7) 


The  canonical  commutation  relation  is  given  by 


[00,  r),  n(f,  r')]  =  -(27r)2<02)(r-  r') 


2  R1R2 


E 


0ik-(r— r') 


(8) 


where  S^2\r  —  r')  is  the  periodic  delta  function  and  k  = 
(si/i?i,  S2/R1)  is  the  2d  momentum  (sj  £  Z). 

To  specify  the  Hilbert  space,  we  develop  the  mode  ex¬ 
pansion  of  the  bosonic  field  </>.  Due  to  the  compactifi- 
cation  condition  (4),  the  bosonic  field  has  the  following 
expansion: 


\  rVi  iN2  <A)  +  tto  t 
d>(t,  r)  =  - x  H - y  H - . 

y  ’  Ri  R2  2tt^/R[R2 


+  0o«cO,  r),  (9) 


where  V12  £  Z  characterize  the  winding  zero  modes  in 
the  x  and  y  direction,  respectively.  The  Fourier  decom¬ 
position  of  the  oscillator  part  </>osc(i,  r)  is  given  by 


4>osc(t,  r) 


1  y _ - _ 

VR1R2  ^  2x/w[k) 

x  [a(k)e_iajt_ik'r  +  at(k)e*“Jt+ik'r]  .  (10) 


According  to  Eq.  (8),  a(k)  satisfies  the  canonical  commu¬ 
tation  relation 


[a(k),at(k,)j  =  5k>k/,  (11) 

where  w(k)  is  the  dispersion  of  the  free  boson  on  a  Eu¬ 
clidean  three-torus  and  given  by 


(12) 


On  the  other  hand,  the  zero  mode  part  satisfies 

[0o,7To]  =  2tt2*.  (13) 


Owing  to  the  27rr  periodicity  of  </>([,  r),  the  eigenvalues  of 
7Tq  needs  to  be  quantized  according  to 


TTO 


ttTVq 

T:\fR\R.2 


Nq  £  Z. 


(14) 


To  summarize,  the  boson  field  <f>(t,  r)  can  be  mode- 
expanded  as 


00, r) 


0o  Apt  ^  1N1  ^  ^  r7V2?) 

2t r^Rjh  2ri?ii?2  R1X  R2  V 

1  1 
+  VR1R2  ^  2v/)J(k) 

x  [a(k)e-iut~ik  r  +  a\k)eiut+ik-r]  .  (15) 
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The  Hilbert  space  Ho  consists  of,  for  each  winding  sector 
specified  by  Ni>2,  the  zero  mode  part  and  the  bosonic 
Fock  space  for  each  k  /  0.  States  in  the  zero  mode 
part  are  labeled  by  the  eigenvalues  of  ttq,  and  hence  by 
iVo.  Furthermore,  different  winding  sectors  are  summed 
over.  In  the  following,  the  part  of  the  partition  function 
associated  to  the  summation  over  iVo,i,2  is  called  the  zero 
mode  sector. 

A.  Modular  transformations  on  T3 

We  now  consider  the  theory  put  on  a  generic  flat  torus. 
A  flat  three-torus  is  parameterized  by  six  real  parameters, 
-Ro, 1,2  and  a,0,j.  For  a  flat  three-torus  T3,  the  dreibein 
is  given  by15 

/  Ro  0  0  W  1  0  0  \ 

eAll  =  0  iZr  0  -a  1  0  (16) 

V  0  0  R2  J  \  -7  -0  1  ) 

where  Rq,  Ri,  and  i?2  are  the  radii  for  the  directions  r, 
x ,  and  y,  and  a,  0,  and  7  describe  the  angles  between 
directions  r  and  x,  x  and  y,  and  r  and  y,  respectively. 
The  Euclidean  metric  is  then  given  by 

A  B  £ 

^  i/OAB 

Rq  +  a2R2  +  7 2R|  — aR\  +  /3y R2  — jR%  \ 
-aR2  +  0  7R2  R2  +  02R2  -0R2  , 

-7 R\  -0R\  '  Rl  ) 

(17) 

The  group  SL( 3,Z)  is  generated  by  two  transforma¬ 
tions: 

/  0  0  1  \  /  1  1  0  \ 
t/r  =  1  0  0  ,  U2  =  0  1  0  .  (18) 

\  0  1  0  /  \0  0  1/ 

Under  the  U2  transformation,  the  metric  is  trans¬ 

formed  as 

9^^  (U2gUZ)^  (19) 

which  corresponds  to  the  changes 

a—>a  —  l,  7— >"f  +  0,  (20) 

while  Rq,  Ri,  R2,  and  0  are  unchanged. 


On  the  other  hand,  Ui  can  be  decomposed  as 

/  0  -1  o\  /  1  0  0  \ 

Ui  =  U[M,  U[  =  1  0  0  M  =  0  0  -1 

\  0  0  1  /  \  0 1  0  y 

(21) 

where  U[  corresponds  to  the  90°  rotation  in  the  r  —  x 
plane  and  M  is  the  90°  rotation  in  the  x  —  y  plane.  The 
generator  U[  acts  on  the  metric  as 

9^^  (U[gU[T)^  (22) 

which  corresponds  to  the  changes 

i?o  ->  Ro/\t\,  Ri  -A  Ri\t\,  n  ->  -ri/|r|2, 

7  -A  —0,  0  — >  7  (while  R2  is  unchanged),  (23) 

where  we  have  introduced 

r  =  a  +  ir0 1,  r0i  =  Ro/Ri-  (24) 

Observe  also  that  under  R0  -A  -Ro/M  and  R\  —A  i?i|r|, 
t2  —A  t2/\t\2.  Hence,  U[  induces  r  —A  —  1/r. 

The  two  transformations  U[  and  U2  correspond  re¬ 
spectively  to  modular  S  and  T^1  transformations  in  the 
t—x  plane,  generating  the  SL( 2,  Z)  subgroup  of  SL( 3,  Z) 
group.  Combining  with  M,  they  generate  the  whole 
SL( 3,  Z)  group.  In  the  following,  we  call  U[M  as  S  trans¬ 
formation  and  U2  as  T_1  transformation. 

Moreover,  the  two  generators  of  SL{ 3,Z),  the  XJ\  and 
U2  transformations  defined  in  Eq.  (18),  satisfy36,44 

IhUl  =  Uf  =  R6  =  (U1R)4  =  (RU1)4  =  1,  (25) 

(! UiR 2)4  =  (R2Ui)4  =  (U1R3)3  =  (R3Ui)3  =  1,  (26) 

(UiR2Ui)2R2  =  R2(UiR2Ui)2  mod  3,  (27) 

where 

R  =  (U2U1)2U^1U2U^1U1U2U1.  (28) 


B.  The  partition  function  on  T3 

In  this  section,  we  calculate  the  partition  function  of 
the  compactified  free  boson  theory  on  the  three-torus  in 
the  presence  of  the  generic  flat  metric.  The  Euclidean 
action  is  given  by 
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S  = 


1 


2-7T 


(2tt)2 

1 

(2^F 


d^O^glg^d^d^ 


2ttR0  2ttR1  2ttR2 


d,T 


dx 


dy 


{dAY  + 


2  /  a2!?2 


i?2 


^2 — b ij  {9x4>Y  +  (a/3  +  T')2  +  + 1 )  (dyY)2 


+2ajr  (dr4>)  {9X<I>)  +  2^  (a/3  +  7)  (dT</>)  (<9y<^)  +  2  (  ^^a(af3  +  7)  +  Yr f3  )  (dx4>)  (dy4>) 

Kq  ti0 


/  R1R2 

\~W 


R2 

Ri 


(29) 


where  0  <  <  27t  are  angular  variables  and  we  noted 

\/\g\  =  RqRiR2,  t  =  Rq9°,  x  =  i?!#1  and  y  =  f?2#2.45 
The  inverse  metric  g Rv  =  (e^)T( 5AB{e*g)  (e\  is  the  in¬ 
verse  of  eA )  is  given  by 


/ 


<r  = 


1 

Rl 

Oi 

a. 

Rl 

“  4- 

1 

a/3+7 

R* 

a(o:/3+7)  , 

Rl 

Rl  + 

Ri 

R2  ' 

(a/3+7)^  ,  /32 

a/S+7 

a(a/3+7) 

,  P 

HI 

*8 

+  R ? 

d2  1  d2 

K0  Jrt1 

\ 


JL 

Ri 


1 

M  ) 


(30) 

In  the  operator  formalism,  the  partition  function  cor¬ 
responding  to  the  action  (29)  is  given  by  the  trace  of  the 
thermal  density  matrix  exp(—2nRoH')  over  the  Hilbert 
space  Rq: 


(where  i,j  are  not  summed  in  RiRjg IJ).  I.e., 


H'  = 


dxdy 

w? 

aR\ 


W 


(dxcj)  +  P^dyf/})2  +  ( dv4 i)2 

'll 


Ro 


+2i:Y5-Tldx(j)  +  2i(af3  +  7)  ~  ndyV 

tin  tin 


(35) 
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The  mode  expansion  for  the  bosonic  field  <fi  is  still  given 
by  Eq.  (15),  where  the  energy  spectrum  w(k)  is  now  given 

by 


Z  =  Tr 


'Ho 


0—2'kRqH' 


(31) 


The  “boosted”  Hamiltonian  H'  and  the  (untwisted) 
Hilbert  space  Ho  are  specified  as  follows.  The  boosted 
Hamiltonian  H'  consists  of  the  “unboosted”  Hamiltonian 
H  (with  a  =  7  =  0),  and  the  momentum  Px,y,  which  in¬ 
duces  the  boost  in  x  and  y  directions,  respectively: 


H’  =  H  +  i^Px  +  i(a/3  +  7)  ^Py 


Rn 


Ro 


(32) 


u/(k)  =  yj  G^kikj 


=  \  9tJslSj  = 


■  (36) 


where 


a  =  j  ^[iP  +  cridiMt], 
Pi  =  J  2n dify,  i  =  x,y. 


and  GlI  is  defined  as 

Y11  gi2 


9 ‘3  = 


921  9  22 



The  Hilbert  space  Ho  is  given  as  a  direct  product  of  the 
(33)  bosonic  Fock  spaces  each  built  out  of  a  given  zero  mode 
state  specified  by  JVo,i,2- 


Next,  we  proceed  to  compute  the  partition  function 
and  study  its  properties  under  modular  transformations 
of  the  three-torus.  The  partition  function  can  be  split 
into  the  zero  mode  part,  which  we  call  Zq,  and  the  oscilla- 
(34)  tor  part,  which  we  call  Zosc.  The  total  partition  function 
is  Z  —  ZqZosc. 
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The  partition  function  of  the  zero  mode  part  is 


Z0=  E  exp 

No,1,2(E^ 


nT2  Ar2  o_„2  D_.fAT  flAM 2  27rr  RqRin2 


Nq  -  2tty2R2t2(N1  +  !3N2y  - 


2t2R2 

2niTiNo(Ni  +  /3N2 )  +  2m')N()N2 


Ri 


where  we  recall  t2  =  Rq/R\,  t\  =  a  and  r  =  t±  +  ir2. 

On  the  other  hand,  for  the  oscillator  part,  the  Hamiltonian  is 


Ksc  =  E 

Mo 


.  Ri 


.Re 


w(k)  +  ia—ki  +  i(a/3  +  7)— fc2 

Uq  -Ko 


a' (k)a(k)  +  E0, 


where  E0  is  the  ground  state  energy  and  needs  to  be  properly  regularized: 

1  r~  \/ det(o)  \  -  1 


e„=  y. 


?•>«,  =  V 


sez2/(o,o) 


2  ^  I  g’isisA2 

sez2/(o,o)  |y  J 1 


l  _  e-2iri?0(^(k )+ia^) 

-1 

1 

L 

tj(t) 

The  partition  function  of  the  oscillator  part  can  be  decomposed  into  the  product  of  the  partition  functions 
dimensional  non-compact  bosons  with  “mass”  given  by  s2.  When  the  “mass”  s2  =  0, 


ZS2=o  =  e-2nRoEo{s2=°) 
where  r/(r)  is  the  Dedekind  eta  function  and  r  is  the  2-dimensional  modular  parameter: 

OO 

V(r):=e^J[(l-qn),  q:=e2™. 

n—1 

On  the  other  hand,  the  other  massive  part  equals  to 

Zs^0  =  e-2vRoEo(»2^0)  1  _  e-2wflo(w(k)-Ma^-HM+7)||) 

=  n  er^,7S2i(r’§s2)’ 

s2ez+  '  2  ' 

where  the  massive  theta  function  @[0S2i7S2](t,  ^ s2 )  is  defined  as 

elaM(r,m)  =  e4^T2A(m^  n|  ^  g— 2'KT2-\i/m2 -\-(n-\-a)2 -\-2niri  ( n-\-a)-\-2nib 

n£Z 

where 

1  _  i  r°° 

A (m,  a)  =  -  ^2  \J m2  +  (re  +  a)2  —  -  /  dk(m2  +  k2)1^2 

2  nez  2 

Thus  the  partition  function  for  the  oscillator  part  equals  to 

1  2 

^ osc  ZS2~  qZ32^o 


1 

V  (t) 


n 0 


s2ez+ 


Ri 

[/3s2,7s2]  l  '  ’  R2 


Together  with  (37),  we  have  completed  the  calculation  of  the  total  partition  function,  Z  =  ZqZosc. 


(37) 


(38) 


(39) 
of  one- 

(40) 


(41) 


(42) 


(43) 


(44) 


(45) 


It  is  instructive  to  compare  the  above  partition  func¬ 
tion  with  the  partition  function  of  the  (l+l)d  compact- 


ified  free  boson.  Performing  dimensional  reduction,  by 
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taking  R2  =  1,  N2  =  0  and  S2  =  0,  the  partition  function 
reduces  to 


Z  = 


1 

Mr)  I2 


E  exp( 

N0,iez 

f  2t:itiNgNi  ) . 


^Nl-2^\2Nl 


(46) 


It  is  also  straightforward  to  show  that  the  oscillator  part 
is  invariant  under  M  transformation  and  thus  the  total 
partition  function  is  invariant  under  M  transformation. 

Finally,  it  is  also  easy  to  check  that  the  partition  func¬ 
tion  is  invariant  under  U2  transformation.  Hence  the 
partition  function  is  invariant  under  the  SX(3,Z)  trans¬ 
formation. 


This  is  the  partition  function  for  the  compactified  free 

boson  in  (1  +  l)d.  III.  THE  SURFACE  THEORY  OF  THE  (3+l)D 

BF  THEORY 


1.  Modular  invariance 


A.  Bulk  and  surface  theories 


We  now  show  that  the  total  partition  function  is  in¬ 
variant  under  the  SL( 3,Z)  transformations. 

For  Z0,  under  U[  transformation,  by  using  the  Poisson 
resummation  formula  twice,  we  have 

Zq{t)  Z0(-1/t)  =  \t\Z0(t)  (47) 

where  the  Poisson  resummation  formula  is 

=  _L  (48) 

n£Z  * 

For  Zosc  part,  under  U[  transformation,  the  massless 
S2  =  0  component  will  contribute  a  1/|t|  prefactor.  The 
massive  part  is  invariant  under  U[  transformation,  since 
the  massive  theta  function  satisfies 


a.  The  bulk  field  theory  The  (3+l)-dimensional  one 
component  BF  theory  is  described  by  the  action 


& bulk  — 


/  M 


K 

— b  A  da  —  a  A  Jqp  —  b  A  Jqv 

Z7T 


J  d^x  dxap  -  apj^p  - 


(52) 


where  n,  v, . . .  =  0, . . . ,  3,  a  =  a^dx11  and  b  = 
(1/2 )bp,vdx^’dx1'  are  one  and  two  form  gauge  fields;  M 
is  the  bulk  spacetime  manifold.  The  “level”  K  is  an  inte¬ 
ger.  The  three  form  Jqp  and  two  form  Jqv  represent  cur¬ 
rents  of  zero-dimensional  (point-like)  quasi-particles  and 
one-dimensional  quasi-vortex  lines,  respectively.  The  BF 
theory  furnishes  the  following  (bulk)  equations  of  motion 


=  @[6  -a]  (r,  (49) 

Thus  the  total  partition  function  is  invariant  under  U[ 
transformation. 

Under  the  M  transformation  which  is  basically  a  7t/2 
rotation  in  the  x  —  y  plane,  the  partition  function  for  the 
zero  mode  part  becomes 


M 


Zn  = 


E 


exp 


Nn 


^0 

2r2R1R2  0 


— 2tti2  RqRiR2 


—  2niaRo[f3Ni  —  N2]Nq  —  2mr)N\NoSJ  . 

(50) 


do  Jqv ,  0  db  —  Jqp-  (53) 

Z7T  Z7T 

The  BF  theory  implements  a  non-trivial  fractional 
statistics  between  quasiparticles  and  quasivortices.  To 
see  this,  we  consider  the  following  configuration  of  quasi¬ 
particles  and  quasivortices: 

Jqp  =  83(C),  Jqv  =  52(S).  (54) 

Here  C  and  S  represent  the  one-dimensional  wold-line  and 
the  two-dimensional  world-sheet  of  quasiparticles  and 
quasivortices,  respectively;  SD-n(JJ)  is  the  delta  func¬ 
tion  ( D  —  n)-form  associated  a  submanifold  J\f  C  M, 
where  D  —  n  =  dimAI  —  dim  Af.  By  definition,  for  any 
n-forrn  An, 


SD—n(Af)  A  An 


Hence,  for  example, 


(55) 


Therefore  the  invariance  of  the  zero  mode  part  of  the 
partition  function  can  be  seen  from  relabeling, 


Jqp  An  — 


A  b  = 


(56) 


f  N\  (  N2 

M  :  l  N2  -A  <  Nt 

[a0  [n0 


Some  useful  properties  of  the  delta  function  forms  are 
summarized  in  Appendix  B. 

In  the  presence  of  these  quasiparticles  and  quasivor¬ 
tices,  we  now  integrate  over  a  and  b  to  derive  the  effective 


action  for  Jqp  and  Jqv.  Since  the  theory  is  quadratic,  this 
can  be  done  by  solving  the  equations  of  motion.  These 
equations,  up  to  a  closed  form,  are  solved  by 

2tt  27T 

b  —  yrrd  Jqpi  ®  Jqv-  (57) 

IV  IV 

(If  the  spacetime  is  trivial,  by  the  Poincare  lemma,  a 
closed  form  is  exact.  If  so,  such  exact  term  does  not  affect 
our  final  result  since,  for  an  arbitrary  closed  submanifold 
J\f,  f  S(Af)(dcj))  ~  f  S(dAf)(j)  =  0.)  From  the  formula 
(B4),  d~1Jqp  and  d~1Jqv  are  determined  as 


b.  The  boundary  theory  On  a  closed  manifold,  the 
BF  theory  is  invariant  under  gauge  transformations  a  — >■ 
a  +  dip,  where  ip  is  zero  form,  and  b  — >■  b  +  where  (  is 
one  form.  In  the  presence  of  a  boundary  (surface),  there 
may  appear  gapless  degrees  of  freedom  localized  on  the 
surface.  The  action  describing  the  boundary  degrees  of 
freedom  can  be  inferred  by  adopting  the  temporal  gauge 
a0  =  ho  =  bot  =  0  (i  =  1,2,3),  solving  the  Gauss  law 
constraints  e^kodkbij  =  e°^kdjak  =  0  by  aj,  =  dkT, 
bij  =  <9iCj  —  djQ,  and  then  plugging  these  back  to  the 
action.  The  resulting  action  is46-48 


d  1  Jqp  =  S2(T>),  where  8fD  =  C, 
d~1Jqv  =  <b(V),  where  dV=S.  (58) 


where  the  two-dimensional  manifold  T>  and  the  three- 
dimensional  manifold  V  are  chosen  such  that  dD  =  C 
and  dV  =  S.  They  are  not  unique,  but  different  choices 
lead  to  different  d~lJqp^qv  which  differ  by  closed  forms. 

Substituting  these  solution  into  the  action, 


2tt  f 

Sbulk  =  /  (d  J qv )  A  Jqp 


=  -^Lk(S,C), 


(59) 


where  Lk  is  the  linking  number  between  Jqp  and  Jqv. 
Hence, 


J  V[a,b}eiSb"lk  =e-^msfi) 

=  qiXjLktSifij)'  (gQ) 

In  the  last  line,  we  assume  the  world-line  C  consists  of  tra¬ 
jectories  Ci  of  many  quasiparticles  each  carrying  charge 
9*  S  Z:  Jqp  =  S3(C)  =  Similarly,  the  world¬ 

line  S  consists  of  trajectories  Si  of  many  quasivortices 
each  carrying  charge  A*  G  Z:  Jqv  =  S2(S)  =  )T7  A i52(Si). 
The  fractional  phase  (when  |K|  >  1)  in  Eq.  (60)  repre¬ 
sents  statistical  interactions  between  quasiparticles  and 
quasi  vortices. 

Once  the  coupling  of  the  gauge  fields  to  the  currents 
is  prescribed,  it  also  specifies  the  set  of  Wilson  loops  and 
Wilson  surfaces  included  in  theory  (see  Eq.  (56)).  If  the 
theory  is  canonically  quantized  on  M  =  K  x  E,  the  set 
of  the  Wilson  loop  and  Wilson  surface  operators  of  our 
interest  is 


exp  im 


exp  in 


(61) 


where  to,  n  are  integers,  and  L  and  S  are  arbitrary  closed 
loops  and  surfaces  in  £.  These  operators  satisfy  the  com¬ 
mutation  relations, 


2m 


fLaJsb  =  ~^I(L’S) 


e imIL  aein  fs  1 


(62) 


where  I(L,  S)  is  the  intersection  number  of  the  loop  L 
and  the  surface  S. 


SdM 


dtdxdy 


K 

2-7T 


tijdiCjdtip  -  V(<pX) 


(63) 


where  i,j  =  1,2.  Here  we  have  added  the  potential 
V(ip,  C),  which  originates  from  microscopic  details  of  the 
boundary  and  is  non-universal.  This  boundary  action 
can  be  obtained  from  the  the  free  scalar  and  the  U{1) 
Maxwell  theories  by  imposing  a  self-dual  (or  an  anti-self- 
dual)  constraint,  e^xdu(^\  =  ±d^ip.47 

The  appearance  of  the  gapless  degrees  of  freedom  on 
the  surface  deserves  more  comments.  In  particular  they 
should  be  contrasted  with  the  gapless  edge  theory  of  the 
(2+l)-dimensional  Chern-Simons  theory.  For  the  single¬ 
component  Chern-Simons  theory  in  (2  +  1)  dimensions, 
the  boundary  is  described  by  the  single-component  chiral 
boson  theory,  which  is  stable  and  cannot  be  gapped  out. 
The  appearance  of  the  gapless  edge  theory  is  necessary 
since  the  bulk  theory  is  anomalous  and  the  anomaly  must 
be  compensated  by  the  degrees  of  freedom  living  on  the 
edge. 

On  the  other  hand,  the  surface  theory  of  the  BF  the¬ 
ory  in  (3  +  1)  dimensions  (also  in  (2  +  1)  dimensions)  can 
be  gapped  out  by  adding  suitable  perturbations  (if  we 
do  not  require  any  symmetry).  In  other  words,  there  is 
no  anomaly  protecting  the  gapless  nature  of  the  surface 
theory.  Nevertheless,  the  appearance  of  the  surface  the¬ 
ory  (63)  can  still  be  understood  in  terms  of  an  anomaly. 
While  the  BF  theory  (both  in  (2+1)  and  (3+1)  dimen¬ 
sions)  is  equivalent  to  the  topological  Zk  gauge  theory, 
which  does  not  have  gauge  anomaly  on  a  manifold  with 
boundary,  the  continuum  action  of  the  BF  theory  ar¬ 
tificially  preserves  U{  1)  symmetry.  The  BF  theory  is 
anomalous  under  t/(l)  in  the  presence  of  a  boundary, 
and  this  anomaly  must  be  canceled  by  gapless  degrees 
of  freedom  living  on  the  boundary.  As  in  the  bulk  the¬ 
ory,  the  boundary  theory  is  invariant  under  the  artificial 
U(  1)  symmetry;  this  symmetry  translates  the  boson  field 
p  — >  p  +  const.  If  the  U(  1)  symmetry  is  strictly  pre¬ 
served,  the  gapless  boundary  theory  cannot  be  gapped, 
as  can  be  inferred  easily  from  the  fact  that  any  gapping 
term  of  cosine  type  cos(mp  +  a)  violates  the  U (1)  symme¬ 
try.  (If  we  use  the  dual  picture  of  the  compactified  boson, 
i.e. ,  the  compact  17(1)  gauge  theory,  the  17(1)  symmetry 
is  equivalent  to  prohibiting  monopoles.)  On  the  other 
hand,  once  we  relax  the  17(1)  symmetry,  which  is,  from 
our  point  of  view,  an  artificial  symmetry  after  all,  this 
gapless  boundary  theory  can  be  easily  gapped  out  by 
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adding  some  relevant  perturbation  and  is  not  stable  at 
all. 

While  this  gapless  surface  theory  is  not  stable  at  all, 
it  does  encode  topological  data  of  the  bulk,  as  we  will 
demonstrate  later.  Let  us  for  now  discuss,  in  more  de¬ 
tail,  the  connection  between  the  bulk  excitations  and  the 
fields  living  on  the  boundary.  In  the  following  we  choose 
Ai  =  S1  x  S,  where  the  spatial  manifold  S  is  a  solid 
torus,  E  =  D2  x  S1,  and  hence  dM.  =  T3 .  Let  us  first 
consider  a  quasiparticle  current  consisting  of  a  quasipar¬ 
ticle  carrying  no  units  of  charges  (no  €  Z): 

jqP{x)  =  n0  J  dTdXdjTh(i)[x  -  X{t)\  (64) 

where  L  is  the  world-line  of  the  quasiparticle,  and  the  co¬ 
ordinate  AM(r)  represents  the  trajectory  of  the  particle. 
For  the  quasiparticle  at  rest,  A'1,2,3(r)  =  const.  =  A'1,2,3, 

j°p{x)  =  n0d(3)(f  -  A).  (65) 

Integrating  the  equation  of  motion  over  the  total  space, 

J  d3X£0ljkdibjk  =  J  d3xj°p  =  n0.  (66) 


Using  Stokes’  theorem,  fs  db  =  fgs  b  = 
(1/2)  fd^bijei:’d2x,  and  substituting  6^  =  diQ  —  djQ, 
this  reduces  to 


f  d  x  —  no- 

Jd  E  K 


(67) 


Hence  adding  a  quasiparticle  in  the  bulk  corresponds  to 
introducing  flux  on  the  surface. 

Similarly,  let  us  consider  to  introduce  a  quasivortex 
source: 

■CM  =  » l -  a»], 

(68) 


where  S  is  the  world-surface  of  the  quasivortex,  and  the 
coordinate  A/J(er)  represents  the  trajectory  of  the  par¬ 
ticle  in  spacetime,  and  n  is  an  integer.  Let  us  consider 
a  straight  quasivortex  at  rest,  stretching  along  a  non- 
contractible  cycle  of  the  bulk  solid  torus.  For  conve¬ 
nience,  this  direction  is  taken  as  the  axdirection  (Fig. 
1).  Then,  j°2  =  j°3  =  0  and 

j°ql(x)  =  n26(x 2  -  X2)S(x3  A3).  (69) 

where  X2’3(a)  =  const,  and  we  have  renamed  the  integer 
n  as  n2.  Integrating  the  equation  of  motion  over  space, 

Li  x  ^  f  dydz £0Uj didj  =  n2  x  Li,  (70) 
2-7T  J 

where  i.  j  =  1,2  and  Iq  =  27ri?!  is  the  length  of  the 
quasivortex  stretching  in  the  ^-direction,  and  we  noted 


FIG.  1.  (a)  The  presence  of  a  point-like  quasiparticle  in 

the  bulk  (solid  torus)  induces  a  fractional  flux  on  the  spatial 
boundary  E  (torus),  (b)  The  presence  of  a  quasivortex  line 
in  the  bulk  twists  the  boundary  conditions  of  the  surface  the¬ 
ory.  Here  and  in  (c),  the  bulk  is  presented  as  a  filled  cylinder 
where  the  top  and  the  bottom  of  the  cylinder  are  identified. 
The  shaded  surface  is  a  sheet  of  the  branch  cut  which  em¬ 
anates  from  the  quasivortex,  and  intersects  with  the  spatial 
boundary  (depicted  by  a  wavy  line).  The  surface  excitations 
experience  a  twisted  boundary  condition  as  they  go  through 
the  branch  cut.  (c)  Similar  to  (b),  a  bulk  quasivortex,  which 
creates  a  branch  cut  on  the  surface  which  now  goes  along  a 
different  cycle  of  the  surface,  is  depicted. 


the  flux  e01y  di  CLj  is  independent  of  x\.  Using  Stokes’ 
theorem,  and  substituting  at  =  diip, 


j  dydyip  =  —n2.  (71) 

Hence  introducing  a  quasivortex  (quasivortices)  along  the 
non-contractible  loop  in  the  bulk  corresponds  to  intro¬ 
ducing  winding  of  the  scalar  boson  on  the  surface. 

One  may  wish  to  develop  a  similar  argument  for  a 
quasivortex  (quasivortices)  stretching  in  the  y-direction. 
(Fig.  1  (c)).  It  should  however  be  noted  that  once  we 
fix  our  geometry  as  above  (Fig.  1  (b)),  loops  running 
in  the  y-direction  are  contractible  in  the  bulk.  In  other 
words,  if  one  constructs  a  solid  torus  by  filling  “inside”  of 
a  two-dimensional  torus,  one  needs  to  specify  one  of  non- 
contractible  cycles  on  the  two-dimensional  torus,  such 
that  after  filling,  this  cycle  now  is  contractible  in  the 
sold  torus. 


B.  The  surface  theory  and  quantization 

We  now  proceed  to  the  canonical  quantization  of  the 
surface  theory.  We  start  from  the  surface  Lagrangian 
density 


C 


~  -  ^^dupdjip. 


(72) 


The  boson  field  tp  is  compact  and  satisfy 


ip  =  ip  +  27 r. 


(73) 
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I.e. ,  physical  observables  are  made  of  bosonic  exponents 

exp[irmp(t,r)],  m€Z,  (74) 

and  the  derivative  of  the  boson  fields  (current  operators) . 
The  winding  number  of  tp  is  quantized,  in  the  absence  of 
bulk  quasiparticles,  according  to 

dx'ditp  =  2 7rlVj,  iVj  G  Z,  (75) 

where  i  =  1,2  and  i  is  not  summed  on  the  right  hand 
side.  On  the  other  hand,  the  gauge  field  Q  is  compact, 
meaning  that  physical  observables  are  Wilson  loops, 

expim  /  m  G  Z,  (76) 

Jc 

where  C  is  a  closed  loop  on  3E  =  T2 .  The  flux  associated 
to  Q  is  quantized,  in  the  absence  of  bulk  quasiparticles, 
according  to 

dxdy  eijdiQ  =  2nN0  (77) 

where  Nq  is  an  integer.  The  canonical  commutation  re¬ 
lation  is 

O'jrj 

Vp{t,r),£ijdiQ{t,  r')]  =  —  <i(2)(r-  r')  (78) 

In  the  following,  we  fix  Ai  and  A2  according  to 


where  the  eigenvalues  of  fa  and  fa^  describes  the  flux 
(associated  with  the  gauge  field  a  in  the  bulk),  and  the 
winding  of  the  tp  field,  respectively.  The  quantization 
conditions  of  these  variables  will  be  discussed  momen¬ 
tarily.  Reflecting  the  compact  nature  of  the  <p  and  Q 
fields,  the  zero  modes  are  compact  variable  =  aM  +  27r 
(/i  =  0,1,2);  For  cto,  the  compactification  condition 
comes  from  the  fact  that  physical  observables  are  given 
as  bosonic  exponents  (74).  Similarly,  for  aii2,  that  phys¬ 
ical  observables  are  given  in  terms  of  Wilson  loops  (76), 
and  that  these  Wilson  loop  operators  must  be  invariant 
under  large  gauge  transformations  imposes  the  compact¬ 
ification  condition,  aq,2  =  aq,2  +  27t. 

From  the  commutator  [p{t,  r ),eijdiQ(t,  r')],  we  read  off 

[a(k),  at(k')]  =  <5k,k', 

From  the  compactification  condition  ao  =  ao  +  27t,  fa  is 
quantized  according  to 

A)  =  ^,  M0gZ.  (83) 

This  quantization  condition  translates  into 

J  dxdy  ei:j diQ  =  .  (84) 


(27t)2 

k2A!A2 


(79) 


This  choice  is  convenient  since  it  gives  rise  to  the  same  en¬ 
ergy  dispersion  as  the  compactified  free  boson  discussed 
in  the  previous  section. 

To  proceed,  we  consider  the  mode  expansion  of  the 
fields.  The  equations  of  motion  are 


-^CijdidtQ  +  1  GtJdtdjp  =  0, 

Z7T  A2 

— ^eikdidtp  +  ^-eikdi(eijdiCj)  =  0.  (80) 

Ztt  a  1 


The  mode  expansion  consistent  with  the  equations  of  mo¬ 
tion  are 


,  s  .  fax  .  fay 

n  r)  =  “0  +  -5-  + 

X £1  1X2 


1 


1 


■£: 


GRAY  2K2A,^;Gk)1/2 


x  [a(k)< 


—  zk- 1 


t(k)-+zkri 


Cj(0  = 


- TV-  4 - „  „  - 

27t  i?7  27ri?ii?2 


-xS. 


j,2 


1 


fa  E 


VR1R2  V  8tt2  ^  u( k)3/2  3 


x  [a(k)< 


—  zk- 1 


t(k)-+zkri 


(81) 


Compared  with  the  quantization  condition  (77),  the  flux 
is  now  quantized  in  the  fractional  unit.  We  will  separate 
Mo  into  its  non-fractional  and  fractional  parts  as 

Mq  =  KiV0  +  no,  N0  G  Z,  n0  =  0, . . . ,  K  —  1,  (85) 
and  write  the  quantization  condition  of  fa  as 

fa  =  No  +  no/K.  (86) 

The  quantization  condition  of  fa<2  can  be  discussed 
similarly.  From  the  commutator  [<p(t,  r),  eijdfaj(t,  r')],  we 
infer 

O'TT'j 

[tijdMt,  r),< fat,  r')]  =  — ^(2)(r  “  O,  (87) 
which  implies 

O7 T1  1 

\fa,a2]-[fa,a1]  =  -  —  (88) 

One  can  choose,  for  example, 

‘Ztt'/  1 

[fa,a2}  =  0,  \fa,ai)  =  —  —  .  (89) 

This  choice  may  be  consistent  with  the  previous  consid¬ 
eration  from  the  bulk  point  of  view,  and  in  particular 
with  the  comment  below  (71).  I.e.,  this  choice  may  cor¬ 
respond  to  choosing  which  non-contractible  loops  on  the 
surface  are  contractible  in  the  bulk,  when  forming  a  solid 
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torus  starting  from  the  two-dimensional  torus  by  filling 
its  “inside”. 

From  the  compactness  of  the  gauge  field  Q,  the  zero 
modes  satisfy  cq  =  a.i  +  27 r,  which  imposes  the  quantiza¬ 
tion  condition 

Mo 

fa  =  -j tz-,  M2  €  Z.  (90) 

Iv 

As  before,  we  split  M2  into  the  fractional  and  non- 
fractional  parts, 

M2  =  K N2  +n2,  N2  €  Z,  n2=  0, . . . ,  K  —  1.  (91) 

With  this,  the  boson  field  obeys  the  twisted  boundary 
condition 

<f(t,  x,y  +  2nR2)  =  ip(t,  x,  y)  +  27r  (E2  +  .  (92) 

While  above  consideration  allows  winding  in  the  y- 
direction  but  not  in  the  ^-direction,  in  computing  the 
partition  functions  of  the  surface  theory  in  the  next  sec¬ 
tion,  we  consider  winding  in  both  directions, 

ip(t,  x  +  27i -Ri,  y)  =  x,  y)  +  2tt  (E  +  , 

<f{t,  x,y  +  2nR2)  =  <p(t,  x,  y)  +  2n  (n2  +  ^  ,  (93) 

That  is 


where  d>  is  a  single  component  boson  theory  compactified 
as  <1>  =  <f>  +  27t,  and  obeys  the  canonical  commutation 
relation  [<&(ar),  ^^^(a:')]  =  (2ni/K)5(x  —  x').  The  zero 
mode  part  of  <f>,  defined  by  the  mode  expansion 

*(*,  X)  =  d>0  -  p(t  +  x)  +  bne~in <*+*> ,  (97) 

n^O 

satisfies  [4?0,p]  =  i/ K.  This  then  suggests  the  quantiza¬ 
tion  rule,  p  =  (integer) /K,  and  the  boundary  condition 
of  the  chiral  boson  field 

$(t,  z  +  2tt)  =  $(t,  x)  —  ^T(il^e§er)  (98) 

Thus,  the  canonical  quantization  naturally  leads  to  the 
twisted  boundary  condition  of  the  chiral  boson  field. 

Quantization  of  the  surface  theory  with  the  above 
twisted  boundary  conditions  gives  the  spectrum  of  lo¬ 
cal  as  well  as  nonlocal  (quasiparticles)  excitations,  which 
obey  untwisted  and  twisted  boundary  conditions,  respec¬ 
tively.  Once  we  specify  the  boundary  condition  (with 
some  integer  vector  n^- o, 1,2)1  the  theory  is  quantized 
within  one  sector  (labeled  by  the  equivalence  class  [n] 
with  the  relation  n  =  n  +  KA  where  A  is  a  vector  with 
integer  entries)  of  the  original  spectrum.  For  this  surface 
theory,  there  are  K3  sectors  in  this  compactified  theory 
and  is  consistent  with  the  K3  ground  states  of  single  com¬ 
ponent  BF  theory  defined  on  T3 . 


fa=i,2  —  Ay2  +  ^1,2/K. 


(94) 


To  summarize,  in  the  presence  of  twisted  boundary 
conditions,  the  mode  expansion  of  the  fields  are  given  by 


r)  = 


Ni 


Hi  K 


Ri 


— 1 


\/RiR2  V  2K2Ai  ^  a/w( k) 


E 


h 


x  [-a(k)e-“(k)‘-ik'r  +  at(k)e+iw(k)t+ik'r 


C.  The  partition  functions  and  modular 
transformations 

Now  we  compute  the  partition  function  (coupled  to  the 
T3  metric): 


Zn°n  i"!  =  Tr« 

rLnQn1n2 


e—2nR0H' 


(99) 


where  7frao„in2  is  the  Hilbert  space  twisted  by  no,ni,n2 
fractional  quantum  numbers,  and 


djdiC j(t,  r) 


N0  +  n0/K 
2itRiR2 


—i 

\J Ri  R‘2 


E 

k^O 


X 


a(k)e-iw(k)t"?;kr 


at(k)e+iw(k)t+ik'r 


(95) 


The  above  consideration  is  somewhat  analogous  to  the 
quantization  of  the  chiral  boson  theory  that  appears  at 
the  edge  of  the  (2+1  )d  Chern-Simons  theory  at  level  K. 
The  (1+1  )d  chiral  boson  theory  defined  on  a  spatial  circle 
of  radius  27t  is  described  by  the  Lagrangian  density 


H'  =  H 


,aRi 
l  —  JTX 

J)  x 

itO 

f  k2A 
H=  dxdv 


i(a/3  +  7)^ 

Ko 

47T2  (  ^  y.  S 

K2Af 


f  K 

Pi  =  dxdy—(eimdiCm)(diy>)- 


Py, 

)2  +  Gij dupdw  , 


(100) 


The  calculation  of  the  partition  function  goes  in  paral¬ 
lel  with  the  calculation  presented  in  the  previous  section 
for  the  free  boson  theory.  To  see  this,  we  note,  from  the 
equation  of  motion, 


eijfaCj 


KAi 

27T 


dt<P 


(101) 


C  = 


K 

47T 


a Mdt 


dx)$, 


(96) 


up  to  a  constant  term.  Thus,  in  terms  of  ip,  the  Hamil¬ 
tonian  density  and  the  commutation  relation  are  given 
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by 

n  =  l^[(dtcpf  +  G^ditpd^\, 

[^t,r),dMt,r'))={^S^\r-r').  (102) 

By  introducing  the  rescaled  field, 

4>  =  K (103) 


the  Hamiltonian  and  the  commutation  relation  can  be 
made  isomorphic  to  those  of  the  free  boson  theory.  The 
compactification  condition  of  the  rescaled  boson  field  is 

r  =  K^E  (104) 

The  partition  function  Zn °"in!  can  now  be  computed 
from  the  partition  function  of  the  free  boson  theory.  The 
zero  mode  part  of  the  partition  function  for  each  excita¬ 
tion  sector  is  obtained  from  Zq  (Eq.  (37))  by  making  re¬ 
placement  Nq  — ►  KA^o+no  and  iVj  — »  Ni+rii) K  ( i  =  1,2): 


g noniri2 


exp{ 

-^0,1,2^^ 


7tK2T2 

2r2i?2 


TV2 


+  2ttitiKNo 


Ni  +  /31V2 


where  we  have  introduced  the  notation 


N »  ■=  N M  +  n^/K.  (106) 

For  the  oscillator  part,  since  the  Hamiltonian  is  the  same 
as  the  oscillator  part  for  the  compact  boson,  the  parti¬ 
tion  function  is  exactly  the  same  as  the  free  boson  case 
presented  above.  Thus  we  have 


-/S2=0-^'S2  7^0 


1 

?7(t) 


n 


s2ez+ 


0 


-l 

[/3s2,7s2] 


(107) 


The  total  partition  function  for  each  sector  is  Zn°nin 2  = 
Zn°n^Zosc. 

Although  the  surface  theory  of  the  (3+l)d  BF  theory 
resembles  the  compactified  free  boson  discussed  in  the 
previous  section,  these  theories  are  physically  different. 


2t~ttR2t2 


N±  +  (3N2 


2  2r27ri?0i?i  ~t2 


R-i 


-Ni 


27 tv 


i7K  N0N2y 


(105) 


For  the  compactified  boson,  the  partition  function  is  in¬ 
variant  under  the  S  and  T  modular  transformations:  It 
is  anomaly-free  and  a  well-defined  theory  on  the  (2  +  l)d 
spacetime  torus.  On  the  other  hand,  for  the  surface  the¬ 
ory,  the  partition  function  for  each  sector  is  not  modular 
invariant  and  thus  it  is  not  a  well-defined  theory  on  the 
(2+l)d  torus.  It  should  be  regarded  as  the  boundary  the¬ 
ory  of  a  higher-dimensional  topological  phase.  There  are 
K3  sectors  determined  by  three  quantum  number  no,  1,2 
and  they  form  a  complete  basis  under  S  and  T  modular 
transformations,  as  we  will  show  now. 

Under  M  transformation,  quantum  numbers  are  trans¬ 
formed  as 


[*1  +  f 

M  :  l  N2  +  %  ->  <  IVi  +  f-  (108) 

[No  +  ^  [iVo+t 

To  discuss  U[  transformation,  we  use  the  Poisson  re- 
summation  to  rewrite  the  summation  over  No  and  N±  in 
Zn°nin2  all(j  rewrite  the  zero-mode  partition  function  as 


gnor 


1 


K| 


E 


exp  • 


_  WT2  m2  -  2r2'7TR2T'2 


N2,Mo,l€% 

27r*riK  Mi 

Irl2  K 


2  r2R 


2T 


Mo 

K 


-  7-/V2 


2nr2R0Ri 

Ri 


Ni 


Mo  Kt  ' 
—  -  jN2 


+  2mpN2Mi  +  2^p-M1 
K 


27rmo 

K 


Mo 


(109) 


Let  us  introduce 


Mi  :=  K N'o  +  n'0t  M0  :=  KiV(  +  n\ 


(110) 
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Then,  the  partition  function  can  be  written  as 
1 


gn0r 


K\t 


Z  Z  exp 

N2,Mo,£Z 

27TZTiK 


7rr2/v2  ~/2  2r27ri?2r2 

2r2i?2|r|2  0  Irl2 
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1V{  -  7lV2  +  2niK(3N2NQ 


N[  -  71V2 

2Trinin'0  2mnon \ 


2  2nr2RoRi 


R‘l 


N2 


K 


K 


(HI) 


where  X),iezK  :=  X)«=o  ■  From  these  expressions,  under 
U[  transformation, 


(JJ[Z)n o’1’"2  =  —  ^  '  e^L{nan'1+n1n'0) 2^n'0n'1n2 

nJ^eZK 


(112) 


Here,  we  made  a  convenient  choice  Ai  =  1/K,  i.e.,  2r2  = 
K.  This  is  the  same  as  the  character  of  the  edge  theory 
of  the  (2+l)d  Zk  gauge  theory  in  its  topological  phase. 
The  effective  Lagrangian  density  of  the  edge  CFT  is  by 

C  =  ^-dt$TKdx$  -  dx$TVdx$,  (115) 


Combined  with  the  M  transformation,  we  can  write  down 
the  modular  S  and  T  matrices: 


S  ,  —  _ A  ,p-T(Bo"2-no’1'i) 

T  /  =  <5  >8  >8  ,  eTpn°ni  11131 

'  nj,n^  uno,nQUni,n1t/ri2,‘n.2C'  • 


This  result  is  consistent  with  previous  works,  Refs.  33  and 
34,  and  also  36,  where  the  action  of  the  modular  trans¬ 
formations  are  calculated  in  the  bulk.  (See  also  other 
related  works:  Refs.  35,  37-39.)  In  terms  of  the  bulk 
physics,  the  S  matrix  describes  the  braiding  phase  be¬ 
tween  particle  and  loop  excitations,  whereas  the  T  ma¬ 
trix  encodes  information  related  to  (3  +  l)d  analogue  of 
topological  spins.33  (See  also  Refs.  35,  37-39.)  The  ex¬ 
act  agreement  between  the  S  and  T  matrices  calculated 
in  the  bulk  and  the  boundary  suggests  there  is  one-to- 
one  correspondence,  the  bulk-boundary  correspondence 
in  (3+l)d. 

The  computed  S  and  T  matrices  (113)  are  expected  to 
be  consistent  with  the  algebraic  relations  in  Eq.  (27): 
As  in  (l+l)d  CFTs,  together  with  the  charge  conju¬ 
gation  matrix  C,  S  and  T  matrices  should  obey  essen¬ 
tially  the  same  algebraic  relations  as  Eq.  (27).  Assum¬ 
ing  the  charge  conjugation  matrix  is  unity,  C  =  1,  we 
have  checked,  for  the  case  of  K  =  2,  3, 4, 5,  the  S  and  T 
matrices  satisfy  all  the  above  constraints  except  the  last 
equation  in  Eq.  (27). 

Before  we  leave  this  section,  as  we  have  done  in  the 
previous  section,  it  is  instructive  to  dimensionally  reduce 
the  partition  functions  of  the  surface  theory  of  the  (3+l)d 
BF  theory.  For  each  given  sector,  after  dimensional  re¬ 
duction,  the  partition  function  is  given  by 


2no,ni 


1 

lrKr)l2 


Z  exp 

No ,  i  E  Z 


-  7TT2K  (iVi  + 

+27rmK  (No  +  (*  +  £)}. 


2 


(114) 


where  K  =  Kcr^  and  V  is  a  symmetric  and  positive  defi¬ 
nite  matrix  that  accounts  for  the  interaction  on  the  edge 
and  is  non- universal.  The  characters  defined  in  Eq.  (114) 
can  be  simplified  as 


Xab(r)  = 



■  \  '  C^5“f°FKi-|-6)2  ^^(Ks+a— K£  — 6)2 


2  /  T 

s.t 


(116) 


where  a  =  no  and  b  =  n\.  There  are  K2  characters  in 
total.  Under  the  S  and  T  modular  transformations,  they 
are  transformed  as 

Xab{r+  1)  =  e2^Xab(r), 

Xab(-l/r)  =  ^  Z  Xa'b'(T)e"2,riiLj^.  (117) 

a',  b' 


D.  Entropy  of  the  boundary  theory 

In  this  section,  we  compute  the  thermal  entropy 

ST  :=  ~^[T  In  X„],  (118) 

obtained  from  the  partition  functions  of  the  boundary 
theory  discussed  above.  Here,  Xa  is  the  partition  function 
in  the  sector  labeled  by  a  =  (no,ni,n2),  and 

1/T  =  2ttR0  (119) 

is  the  inverse  temperature. 

While  St  is  defined  for  a  system  with  a  real  (physi¬ 
cal)  boundary,  it  is  expected  to  carry  information  on  the 
universal  topological  part  of  the  entanglement  entropy 
(the  topological  entanglement  entropy) .  The  latter  is  de¬ 
fined  for  the  bulk  system  (the  BF  theory)  defined  on  a 
manifold  without  a  physical  boundary,  and  obtained  by 
integrating  out  (tracing  over)  a  subregion  B  (compliment 
to,  say,  subregion  A). 8,49-51 
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We  are  interested  in  the  entropy  St  in  the  limit 
Ri/Ro  —>  oo  and  R1/R2  — »  00.  (We  could  also  equiva¬ 
lently  take  the  limit  with  Ri  and  R2  exchanged,  in  which 
case,  we  have  to  resum  differently  but  the  result  would 
be  the  same.)  To  evaluate  the  entropy  in  this  limit,  we 
first  make  use  of  the  5-modular  transformation,52  and 
write 

St=-^[T In  (SbaXb(-l/r))].  (120) 

In  the  above  limit,  only  the  identity  charac¬ 
ter  gives  rise  to  the  dominant  contribution, 
limfli/iJ0^oo,ii1/fl2-).ooXa('r)  =  S°Xo(-l/T),  as  seen 
from  Eq.  (111).  Hence 

St\r1/r0^oo,r1/r2^oo  =  ~j^  [Tin  («S°Xo)]  ■  (121) 

Then  using  the  modular  S  matrix  computed  in  the  pre¬ 
vious  section, 

St  I  Ri  /  Rq — voo  ,Ri/Ft.2 — voo  clT 

=  —  in  K  +  [T  in  x°]  •  (122) 

The  first  term  is  the  subleading  term,  and  identical 
to  the  bulk  topological  entanglement  entropy,  although 
St  and  the  entanglement  entropy  are  defined  differently. 
The  second  term  is  the  extensive  piece,  which  basically 
corresponds  to  the  entropy  of  the  free  boson  and  is  the 
usual  leading  order  term.  (When  St  is  interpreted  as 
the  entanglement  entropy,  the  second  term  corresponds 
to  the  area  law  term.) 

IV.  THE  SURFACE  THEORY  OF  THE  (3  +  1)D 
BF  THEORY  WITH  THE  0  TERM 

Recall  that  in  the  surface  theory  of  the  BF  theory  dis¬ 
cussed  in  the  previous  section,  there  are  three  quantum 
numbers  Mq  -[  2,  which  we  wrote  in  terms  of  their  non- 
fractional  and  fractional  parts  as 

77/ 

MP  =  NP  +  ^,  ^  =  0,1,2.  (123) 

These  quantum  numbers  in  the  surface  theory  can  be 
interpreted  as  arising  from  the  presence  of  bulk  quasi¬ 
particles  or  quasi-vortices;  Mi  .2  represents  the  frac¬ 
tional  winding  number  of  the  ip  field  induced  by  a  bulk 
quaxi- vortex,  whereas  Mo  represents  a  fractionalized  flux 
threading  the  surface  induced  by  a  bulk  quasi-particle. 

In  this  section,  we  consider  the  following  “twist”  of  the 
quantum  number 


in  the  surface  theory  of  the  BF  theory,  where  Qip  are 
fixed  integers,  and  we  have  introduced  the  notation 

QxM  ~Q1M2-Q2Ml.  (125) 

This  twist  can  be  induced  by  considering  a  modification 
of  the  BF  theory  by  introducing  the  0  term  (axion  term) . 
In  the  next  section,  we  will  consider  a  similar  twist  to 
discuss  three-loop  braiding  statistics. 


A.  The  BF  theory  with  the  O-term  in  (3+l)d 


We  motivate  the  twist  (124)  by  considering  the  fol¬ 
lowing  modification  of  the  bulk  BF  theory  by  adding  a 
0-term: 


Sbulk 


—  b  A  da - J-77  dQ  A  a  A  da 

2tt  8tt2 


—  a  A  JqP  —  b  A  Jqv 


(126) 


In  the  second  term  (the-0  term  or  axion  term),  p  is  a 
parameter,  specific  value  of  which  will  be  discussed  later, 
0  is  a  non-dynamical  background  field,  and  we  consider 
an  inhomogeneous  but  time-independent  configuration  of 
0,  which  will  be  specified  later.  Compared  to  the  stan¬ 
dard  form  of  the  0  term,  Qda  A  da,  we  have  done  an 
integration  by  part  and  put  the  derivative  acting  on  0. 
Since  the  0  field  is  non-dynamical,  we  will  interpret  the 
presence  of  the  0  term  as  an  introduction  of  a  static 
defect.  In  Ref.  53,  a  similar  effective  action  has  been 
proposed  to  describe  the  thermal  and  gravitational  re¬ 
sponse  of  topological  defects  in  superconducting  topolog¬ 
ical  insulators.2'  We  also  note  that  the  BF  theory  with 
the  0-term,  Qda  A  da,  has  been  proposed  to  describe 
the  fermonic  and  bosonic  topological  insulators.  In  Ref. 
37,  the  BF  theory  with  the  0  term  was  used  to  discuss 
three-loop  braiding  processes. 

To  see  the  0-term  induces  the  twist  (124),  we  assume 
the  following  configuration  of  the  ©-field: 

Q±x  Q2y  ,  . 

®(x,y,z)  =  (127) 

rti  -K  2 

where  Q\p  are  fixed  integers.  From  the  equation  of  mo¬ 
tion, 


^Xpdxap=j^, 

^vXpdvbXp  +  ~^e^XpdvQdxap  =  jpp.  (128) 

By  plugging  the  first  equation  into  the  second,  these 
equations  of  motion  reduce  to 


M0  -A  M0  + 


Q1M2  —  Q2M1 


K 


=  M0  + 


Q  x  M 


(124) 


K 

27T 

K 

47T 


epvXpdxap=j%, 


epvXpdvbXp  =  - 


P 

2ttK 


K 


dvQjqv  +  Jqp- 


(129) 
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In  the  presence  of  quasiparticle  and  quasivortex  sources, 
(65)  and  (69),  the  equations  of  motion  integrated  over 
space  are 


J  dydz  £01lJ diCij  =  Mi, 

J  dxdz  diaj  =  M2, 

^  J  d3xe0ijkdibjk  =  QiMi  +  N0, 


(130) 


where  Mit2,N0  G  Z  and  we  noted 

dte  J  d3xj°qlv  =  O-Mi  x  (27 TRi)  =  2tt QiMi  (131) 

(*  is  not  summed  over).  These  can  be  reduced  to,  by 
using  Stokes’  theorem, 

J  dydyifi  =  K  N2, 

f  2n 

/  dxdxtp  =  —ATi, 

J  d?xeijdiQ  =  ^-N0  +  ^-(Qx  N),  (132) 

where  Mi  =  — Ara  and  M2  =  N±.  Hence,  upon  choosing 
p  =  1,  in  the  presence  of  the  defect  field  0,  the  quantum 
number  N0  in  the  surface  theory  is  “twisted”  as  in  Eq. 
(124). 

We  observe  that  the  following  action 


SLk  =  J  d*x^e^dvbXpa 


■  N  n  -  i^-\b  - 

JqpaH  Jqv  2  yJ^1' 


(dftdvQ  -  a„dM0) 
Z7Tl\  J 

(133) 


we  note  that  the  quantization  rule  of  the  zero  modes  are 
determined  both  by  (a)  the  canonical  commutation  rela¬ 
tions  and  (b)  the  compactification  conditions.  The  com- 
pactification  condition  is  determined  by  declaring  physi¬ 
cally  observable  Wilson  loop  operators.  This  in  term  is 
determined  from  the  coupling  of  the  theory  to  the  cur¬ 
rent.  Therefore,  in  the  original  theory,  (a)  is  abnormal 
but  (b)  is  normal.  In  the  modified  theory,  (a)  is  normal 
but  (b)  is  abnormal.  In  the  next  section,  we  demonstrate 
this  by  deriving  the  quantization  conditions  (132),  de¬ 
rived  from  the  bulk  point  of  view  here,  in  terms  of  the 
boundary  theory  of  S'bulk.  In  Appendix  A,  we  quantize 
the  boundary  theory  of  the  original  theory,  Sbuik ,  to  de¬ 
rive  the  quantization  rule  (132). 


B.  The  surface  theory  and  partition  functions 


1.  The  compactification  conditions  and  quantization  rules 

We  now  proceed  to  consider  the  surface  theory  of  the 
bulk  theory  (133).  Without  sources,  the  surface  theory  is 
described  by  the  same  Lagrangian  density  as  the  surface 
of  the  BF  theory,  (72),  and  hence  has  the  same  canonical 
commutation  relations.  This  immediately  means  that  the 
oscillator  part  of  the  surface  theory  can  be  treated  in 
exactly  the  same  as  before.  On  the  other  hand,  reflecting 
the  abnormal  coupling  of  the  gauge  fields  to  the  currents 
in  the  bulk  action  S'bulk,  the  compactification  conditions 
of  the  boundary  fields  ip  and  £  are  modified,  as  we  will 
now  discuss. 

As  we  noted  earlier,  the  coupling  to  the  current  can 
be  written,  e.g.,  f  d4x  j^d^  =  §L  a.  Thus,  introducing  a 
proper  current  corresponds  to  introducing  a  Wilson  loop. 
If  we  now  consider  a  Wilson  line  L  that  is  spatial,  and 
that  ends  at  the  boundary, 


shares  the  same  equations  of  motion,  Eq.  (129),  as  the  BF 
theory  with  0  terms,  (126).  Hence,  the  boundary  theory 
derived  from  Sbulk  has  the  same  quantization  rules  of  the 
zero  modes  as  the  boundary  theory  of  Sbuik-  In  the  next 
section,  we  will  consider  the  boundary  theory  derived 
from  S'bulk,  and  its  partition  functions. 

To  contrast  the  two  theories  Sbiuk  and  Sbulk,  we  note, 
in  Sbuik,  that  the  coupling  to  the  currents  are  “normal” 
while  the  commutators  are  “abnormal” ,  in  the  sense  that 
the  commutators  among  fields  a,  b  are  modified  due  to  the 
presence  of  the  theta  term.  On  the  other  hand,  in  Sbulk, 
the  commutators  are  normal  (the  same  as  the  ordinary 
BF  theory)  while  the  coupling  to  the  current  is  “abnor¬ 
mal”.  (Since  that  the  commutators  are  the  same  as  the 
ordinary  BF  theory,  Sbulk  and  the  corresponding  bound¬ 
ary  theory  can  be  analyzed  in  a  complete  parallel  with 
the  BF  theory  -  a  practical  reason  why  we  will  consider 
on  Sbulk  in  the  following  -  expect  for  the  zero  mode  part.) 

In  spite  of  these  differences,  these  theories  lead  to  the 
same  quantization  conditions  (the  same  “lattice”  of  quan¬ 
tum  numbers)  of  zero  modes.  To  see  how  this  is  possible, 


a  =  I  dtp  =  I  p  =  p(dL) 


(134) 


IdL 


where  we  noted  dL  is  a  point,  and  we  have  used  the 
solution  to  the  Gauss  law  constraint,  da  =  da<p  (a  = 
1,2,3).  Thus, 


expim  J  d  =  exp  [imtp(dL)]  (135) 

This  means  that  ip  is  compactified  with  the  radius  27t. 

Let  us  repeat  the  same  exercise  for  the  coupling  to  the 
quasi  vortex  current: 


dAxrv\ 


h  p 


is  L 


qv  2  L 

b  —  d  A  d<d 
2ttK 


{d^dv®  -  dvd^Q) 


(136) 


where  S  is  the  world  surface  of  a  quasivortex  (quasivor¬ 
tices).  In  the  presence  of  a  boundary  and  using  b  =  d(, 
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this  is  evaluated  as 


Inverting  this  relation, 


is  L 


dC-  -^-dtpAdQ 
s  2ttK  ^ 


IdS  L 



(137) 


where  the  boundary  of  the  world  sheet  is  on  the  surface. 
We  thus  have  a  Wilson  line  on  the  surface: 


exp  im 


ids  L 


c- 


2ttK 


if  A  d<d 


(138) 


Pp  =  <iu£  =  (5()  (148) 


where 


u0  =  (1,0,0)J 


*  =  ,^.i,o)  -  ■  (149) 


We  now  consider  the  case  where  dS  is  along  the  x-  or 
y-  cycles.  Recalling  the  mode  expansion  Eq.  (81),  and 
noting  (p/27rK)  fL  fd&  =  (p/K)a0 Qi  the  zero  modes 
enter  into  the  integral  Jgs  [£  —  A  d0]  through  the 

following  combinations 

cki  —  ^Qicio;  a2  —  ^-Q2Q:o-  (139) 

Together  with  ao ,  the  following  three  linear  combina¬ 
tions 


Renaming  the  integers  as  mo  — >  No,  mi  — >  —  IV2,  and 
?ri9  — >  N\,  Eq.  (148)  is  nothing  but  the  quantization  rule 
(132). 

2.  The  partition  functions 

With  the  twist  (124),  we  can  now  write  down  the  zero 
mode  partition  function.  Let  us  recall  the  partition  func¬ 
tion  of  the  BF  surface  without  the  theta  term,  Z"0"' n'2 , 
defined  in  Eq.  (105).  For  later  use,  we  write  Zn 0"i”2  as 


v a  =  0,1,  2  (140) 


are  angular  variables,  where 
v°  =  (1,  0, 0)T, 


Noting  the  commutation  relations  among  zero  modes, 


Znonin2=  /K  {Mo  Mi) 

-No, 1,2^ 

=  Y,  UiKNo  +  no^Ni  +  m),  (150) 

No,  1,2 

where  /k  is  defined  by  the  summand  in  Eq.  (105),  and 
recall  M M  =  KN^+n^.  We  will  call  the  partition  function 
resulting  from  the  twist  Zq^  q1"2 .  It  is  given  by 


[&o,Po]  —  [a  1,  -/? 2]  —  [a2,P\]  —  — ,  (142) 

we  consider  the  linear  combinations 

Kh,  P  =  (Po,~P2,Pi),  a  =  0,1, 2,  (143) 

where  wa  are  translation  vectors  reciprocal  to  va: 

«  =  *6°.  (144) 

Explicitly,  they  are  given  by 

W°  =  (1’  1  Wl  =  (°> l- °)>  W2  =  (°i  °>  !)■ 

(145) 

Then,  in  the  “rotated”  basis,  the  commutation  relation 
takes  the  following  canonical  form: 


ZqTqT=  E  fK(KN0+no  +  ^-^,KNi+nX 

iv0,i,2ez  '  ' 

(151) 

To  proceed,  we  write 

Mi  =  K Ni  T  Tii  =  K2Ni  +  Rti  +  Tii , 

Qi  =  KRi+ri  =  K2Ri  +  Ksi+ri,  (152) 


where  new  integers  Ni,Ri,Ri  and  Zk  variables  ti,  s,;,  rt 
are  introduced.  In  the  following,  we  will  show  that  the 
zero  mode  partition  function  depends  on  Qi  only  through 
Ti,  and  hence  can  be  denoted  as  z5”°r"in2,  and  that  the 
partition  function  can  be  written  as 

Z™°£n2  =  Y  X"-0"1”2,  (153) 


K&P,  WbPu]  =  =  ^<5£.  (146) 

Due  to  the  compacticity  of  wbPv  takes  on  values 

w%fiv  =  —  x  mb,  mh= ij2,3  €  Z.  (147) 


where  we  have  introduced 

nn  :=  no  +  s  x  n  +  r  x  t  mod  K, 

{no  :=  Kn0  +  (r  x  n) 

Hi  :=  K t\  +  ni 

n2  :=  Kf2  +  n2, 


(154) 
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(i.e.,  n0  =  (no  +  s  xn  +  rx  t)% K)  and  Xn°nin2  is  defined 
by 


j£n0r. 


E  /k(: 


A0,Nli2ez 


KA0+^,K2N.l+n. 


(155) 


To  show  Eq.  (153),  we  start  by  writing  the  partition 
function  in  terms  of  variables  introduced  in  Eq.  (152): 


ynoniri2  _ 

Qi=K2  Ri+Ksi+n 


=  E  E 

n0,  JVii2ezti,2eZK 

x  /k  +  no  +  ^  ,  K“iVj  +  Kti  +  n^j  .  (156) 

By  further  introducing 

«o  =  no  +  s  x  n  +  r  x  t, 

Aq  =  Nq  +  (K  R  +  s)xN  +  Rxn  +  rxN,  (157) 
and  noting  the  equality 

KN0  +  „0+^A=KA  <,  +  ao  +  rAA  (158) 


Then, 


ynonin2 

ZjQi=K2Ri+Ksi  + 


K 


,,  =  E  E 

tl,2£^K 


K 


x  /k  ^KA0  +  no  H  jv — ,  K2iVj  +  Ktj  +  n^  )  .  (159) 

We  now  fix  t\  2  and  consider 


■ynQn\ri2  ,tit2 
^Qi=K2  Ri+Ksi+n 


=  E 

No,N12(z2 


x  /k  ^KA0  +  ao  d - — ,  K2iVj  +  Ktj  +  riij  .  (160) 

Note  that  once  £1,2  are  fixed,  ao  is  fixed.  Converting  the 
summation  over  Nq  to  a  summation  over  Aq , 


■y-noniri2  ,tit2 
^  Qi=K2  Ri+Ksi  + 


r=  E 

-Aq  5-^1 ,2 


x  /k  ^KAo  +  ao  d - jV — ,  K “Ni  +  d -  riij  .  (161) 

Note  that  the  s  and  Uq  dependence  of  the  right  hand 
side  comes  only  from  ao.  Also,  after  converting  the  sum 
12  n0  Ea0i  R*  dependence  is  gone.  So,  we  write 
Xq° =K2R*+Ksi  +n  simply  as  Observe  that 

X'q  ’~k2r'+Ks  +r  aPPears  to  depend  on  nine  Zi<-valued 
parameters,  no, 1,2,  £1,2,  si,2,  7T,2  After  the  reorganization 
we  have  just  done,  we  lost  s  i,2,  and  we  now  only  have 
six  Zk  parameters,  7x1,2,  £1,2,  tt>2  and  ao.  While  ao  is  not 
Zx-valued,  we  can  shift  ao  such  that 

(162) 


where  the  second  term  takes  values  0, ...  K  —  1.  Then, 

v-[ao]nin2,tit2  _  \  ' 

rir2  /  y 

A0,Ni,2eZ 

x  /k  ^KA0  +  [ao]  d - — — ,  K2W  -|-  Kti  d-  n-i'j  ■ 

(163) 

Observing  that  A^r"1"2’*1*2  depends  on  [ao], 7x1,2, ti,2 
only  through  710,1,2  defined  in  Eq.  (163),  rewriting  Eq. 
(163)  in  terms  of  no  1  2  completes  the  derivation  of  Eq. 
(153). 


3.  Modular  transformations 

We  now  discuss  the  modular  properties  of  the  partition 
functions.  Under  the  U2  transformation,  the  zero-mode 
partition  functions  are  transformed  according  to 


(C/2  X)n°nin2  =  e~  K^n°ni xn°nin2 


(164) 


On  the  other  hand,  under  the  U[  transformation,  the 
partition  functions  are  transformed  as 


/  Y\nonm2 


(U[X) 


1 

MK2 


E 

n'n  !  ,r' 


g  fjT  (nO n'l  +nl n 0 )  XnO nl ”2 
(165) 


where  n'0  =  Kn.g  +  r(n2  —  t^tx^. 

Observe  that,  upon  the  U[  transformation,  partition 
functions  with  new  parameters  r\  and  r'2  generated. 
Since  r\  and  r2  are  the  given  quantum  numbers  from 
the  O  term,  the  action  of  modular  transformations  is  not 
closed. 


V.  COUPLING  TWO  BF  THEORIES  - 
THREE-LOOP  BRAIDING  STATISTICS 

In  the  twist  (124),  the  integers  Qi,2  are  fixed  and 
treated  as  a  background.  I.e.,  O  is  a  non-dynamical  field. 
We  have  seen  that  the  surface  partition  functions  do  not 
form  a  complete  basis  under  modular  transformations. 
To  circumvent  this  issue,  one  may  consider  to  treat  Q\ 
and  Q2  as  dynamical  variables,  which  may  come  from 
another  copy  of  the  BF  theory.  In  this  section,  we  will 
discuss  two  copies  of  the  BF  surface  theories  which  are 
coupled  via  cubic  terms. 

Let  us  start  from  two  decoupled  copies  of  the  BF  sur¬ 
face  theories.  Let  M0, 1,2  and  Q0,i,2  label  different  twisted 
sectors  of  the  first  and  second  copy,  respectively.  We  con¬ 
sider  to  twist  these  quantum  numbers  by 


M0  -A  M0  + 
Qo  -A  Qo  + 


Q  x  M 
K  ’ 
M  x  Q 


ao  =  K  x  (integer)  +  [ao] 


K 


(166) 
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Here,  unlike  Eq.  (124),  both  M  and  Q  are  dynamical 
variables. 

In  the  next  section,  we  start  by  introducing  an  (3+l)d 
bulk  field  theory,  Eq.  (167),  or  its  alternative  form  (173), 
which  realizes  precisely  the  twist  (166).  We  will  analyze 
the  modular  properties  of  the  resulting  zero  mode  parti¬ 
tion  functions  at  the  surface.  The  oscillator  part  of  the 
partition  function  is  simply  given  by  the  partition  func¬ 
tion  of  the  two  decoupled  copies  of  free  boson  theories. 
By  computing  the  S  and  T  matrices  acting  on  the  zero 
mode  partitions,  we  argue  that  the  action  (167)  realizes 
three  loop  braiding  statistics.  While  we  were  finalizing 
the  draft,  a  preprint  54  appeared  where  the  similar  bulk 
actions  were  discussed  and  conjectured  to  realize  three- 
loop  braiding  statistics. 


A.  The  bulk  field  theory 


a.  The  cubic  theory  Let  us  motivate  the  twist  (166). 
We  propose  to  work  with  the  following  bulk  action: 


Sbulk  — 


K 

27T 


Sijb1  A  daJ 


+  — -a1  A  a2  A  da 2  +  -^a2  A  a1  A  da 1 

47T2  47T2 


—  Sub1  A  Jjqv  —  djja1  A  Jqp  , 


(167) 


where  /,  J  =  1,2  and  pi^  are,  as  the  level  K,  constant 
parameters  of  the  theory.  Similar  action  has  been  dis¬ 
cussed  in  Ref.  55-57.  The  equations  of  motion  are 


K  j  j 

—da  =  Jqv , 


K 


Pi  2 


— — db  +  — — t?  cl  A  da 
2  7T  47T2 

-toa2Adal]toda2Aal  =  jl 
2n2  47T2  qp 

— db 2  +  -^a1  A  da 1 
27T  47T2 

-  A  da2  +  7-^da1  A  a2  =  J2 

27t2  47t2  qp 


(168) 


As  in  our  previous  discussion  in  the  BF  theory  with 
and  without  the  theta  term,  let  us  consider  a  fixed,  static 
quasiparticle  and  quasivortex  configuration  and  integrate 
the  equation  of  motion  over  space.  By  solving  the  first 
equation  of  motion  as  a1  =  (27r/K)(d_1  Jqv),  plugging  the 
solution  to  the  second  and  the  third  equations  of  motion, 


and  integrating  over  space, 

+  1Jqv )  A  Jqv  +  Jqp’ 

+  %  j  (d~ljqV)  A  J2qv  +  J2p,  (169) 

where  note  that  in  the  static  configurations  considered 
here,  Jqv  is  a  delta  function  one  form  supporting  a  spatial 
loop,  whereas  Jqp  is  a  delta  function  three  form  support¬ 
ing  a  spatial  point.  Correspondingly,  d~1Jqv  is  a  delta 
function  0  form  supporting  a  three  dimensional  manifold. 
The  contributions  to  fs  db1  is  coming  from  quasivortex 
loops,  /E(d_1  Jqv)A  Jqv,  are  given  in  terms  of  their  linking 
number. 

Considering  now  the  specific  geometry  E  =  D2  x  S 1 
with  the  boundary  (surface)  dT,  =  T2,  we  can  derive 
the  quantization  rule  of  the  zero  modes  of  the  bound¬ 
ary  fields.  Using  the  Gauss  law  constraint  to  write  the 
boundary  conditions  in  terms  of  if1  and  C)7,  the  bulk  equa¬ 
tions  of  motion  translate  in  to 


[  rfC2  -  -pz  [  dT1  A  dip2  =  Q0.  (170) 

2tt  JdYi  4t r2 

With  Pi  =  P2  =  K,  these  correspond  precisely  to  the 
twist  (166). 

Note  that  if  we  naively  gauge  transform  as  b1  -A  b1  + 
d^1  and  a1  — >  a  +  dip1 ,  we  find  that  the  theory  is  not 
gauge  invariant.  Moreover  these  gauge  transformations 
are  not  generated  by  Gauss  constraints.  We  propose  the 
following  alternative  gauge  transformations: 

b 1  — >  b'1  =  b1  +  d C1  —  — — (a2  A  dip 1  +  dip2  A  a1)  , 
27tK  v  7 

b 2  -A-  b '2  =  b2  +  dC,2  —  (fll  A  dip2  +  dip1  A  a2)  , 
27tK 

a1  -A  a'1  =  a1  +  dip1 .  (171) 

Therefore  the  action  with  cubic  terms  in  Eq.  (167)  is 
gauge  invariant.  On  the  other  hand,  as  for  the  coupling 
to  the  sources,  by  demanding  the  gauge  invariance,  we 
can  read  off  the  conserved  currents,  which  are  modified 
due  to  the  presence  of  the  cubic  terms  and  the  modified 
gauge  transformations. 

On  an  open  manifold,  the  action  picks  up  a  gauge 
anomaly  on  the  boundary  under  these  gauge  transfer- 
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mations 


where 


Rbulk\P  ;  bl  ]  Sbulkfoi  d] 

+  7T  [  hjdC1  A  aJ 
J dM 


+  — -=  /  (pi cbp1  A  dip2  A  a2  +  P2 dip2  A  dip1  A  a1) . 

^2  J  dM 

(172) 


This  anomaly  then  must  be  compensated  by  an  appro¬ 
priate  boundary  field  theory. 

b.  The  alternative  quadratic  theory  Instead  of  tack¬ 
ling  the  cubic  theory  (167)  and  the  corresponding  surface 
theory,  as  in  our  discussion  in  the  BF  theory  with  theta 
term,  we  consider  an  alternative  form  of  the  theory.  We 
note  that  the  equations  of  motion  (168)  can  be  derived 
from  the  following  alternative  action: 


<?'  = 
^ bulk 


2^:  J  dijb1  A  da'J 


P2  1 
.a 


2ttK 


b2  + 

2ttK 


A  a2 
A  a1 


5ijaJ  A  Jqt 

A  Jqv 

A  Jqv. 


(173) 


Unlike  Sbuik ,  this  theory  is  quadratic.  Integrating  over 
a 1  and  b1 ,  one  obtains  the  effective  action  of  the  currents 


J  £>[cf  67]eiS^  =eiS°st  (174) 

where 

sc  ff  =  -  K  jid^J^AJ1^ 

+  Pi  J (d-1  Jqv)  A  (cT1  J2V)  A  J2V 

+  P2  J 0 d ~ljqv)  A  (d  1  Jqu )  A  Jqv.  (175) 

The  first  term  in  the  effective  action  describes,  as  in  the 
ordinary  BF  theory,  the  quasparticle-quasivortex  braid¬ 
ing  statistics  while  the  second  and  third  terms  include 
interactions  among  three  quasivortex  lines. 

From  the  coupling  to  the  currents,  we  read  off  the  Wil¬ 
son  loop  and  Wilson  surface  operators  in  the  theory: 


B1  :=  b1  +  ^d-a1  A  a1 ,  (178) 

27tK 

and  as  before  the  repeated  capital  Roman  indices  are  not 
summed  over.  Note  also  the  triple  commutator  among 
fs  B1  is  computed  as 


Is  B1 1  fs<  BJ\ ,  fs„  BK 

-  4?rpj  (su5Jk  -6rJSJK 


K3 


I(S#S\S").  (179) 


To  make  a  comparison  between  the  cubic  and  quadratic 
theories,  in  the  cubic  theory,  the  canonical  commutation 
relations  differ  from  the  ordinary  BF  theory,  while  they 
remain  the  same  in  the  quadratic  theory.  In  fact,  in  the 
cubic  theory,  the  commutator  among  fields  generates  an¬ 
other  field,  [6,  b]  ~  a,  schematically.  On  the  other  hand, 
the  set  of  Wilson  loop  and  surface  operators  in  the  cubic 
theory  is  conventional  (i.e. ,  identical  to  the  ordinary  BF 
theory)  while  it  is  modified  in  the  quadratic  theory  as  in 
(176).  In  spite  of  these  differences,  the  algebra  of  Wilson 
loop  and  surface  operators  of  the  two  theories  appear  to 
be  identical.  Therefore,  we  argue  that  the  two  theories 
are  equivalent.  In  the  following,  we  will  proceed  with  the 
quadratic  theory. 

c.  the  quantization  rule  of  the  zero  modes  We  now 
derive  the  compactification  condition  of  the  boundary 
fields  from  Eq.(176).  In  the  presence  of  a  boundary  and 
using  b1  =  df1  and  a1  =  dip1 ,  the  surface  operators  re¬ 
duce  to 


exp  im 


PI 

2ttK 


dip1  A  dip1 


=  exp im 


PI 

2ttK 


ip1  A  dip1 


(180) 


where  the  boundary  of  the  world  sheet  is  on  the  surface. 
We  now  consider  the  case  where  dS  is  along  the  x-  or  y- 
cycles  on  the  surface.  Recalling  the  mode  expansion 


</(r)  =  a7 

C/(0  = 


“j 


+ 
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R2 
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2nRn  2ttR1R2 


X§j,  2  +  ■ 


(181) 


exp 

im  /  a1 

,  exp  in 

J  L  . 

Js  L 

A  o1 

2ttK 


,  (176) 


where  n  and  m  are  integers,  L  and  S  are  arbitrary  closed 
loop  and  surfaces,  respectively,  and  we  introduced  the 
notation  1  =  2  and  2  =  1,  and  the  repeated  capital  Ro¬ 
man  indices  are  not  summer  over  here.  These  operators 
(or  rather  their  exponents)  satisfy 


fcaIJsBJ 


2ni 

K 



2ipj 


s  fs#s '  + 


J  i  2*P  J  xiJ 


K2 


fs#S'  aJ 

(177) 


the  zero  modes  enter  into  the  integral  /9S[C7  + 
(pj/27tK)95/  A  dip1]  through  the  combinations 

S7  +  f ' yl  (182) 

We  thus  conclude 



v2aa 2 


(183) 


are  angular  variables,  where 

v10  =  (1, 0,  0)T, 

v11  =  (pi^i/K,  1,  Of,  v12  =  (pj/3|/K,  0,  If.  (184) 
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The  rest  of  the  discussion  is  essentially  identical  to  the 
analysis  made  in  Sec.  IV  B  1.  We  recall  the  commutation 
relations  among  zero  modes 

Wq,Pq]  =  Wn-Pi]  =  [<*2,  £2]  =  jE/j,  (185) 

the  following  linear  combinations  of  the  zero  modes  are 
integer-valued 


Under  the  U2  transformation,  the  product 
^n0n1n2j^r0r1r2  jg  invariant  up  to  a  phase, 

(. u2x)n°nin2{u2xy°fif; 2 

_  g  — ^(nom+rori)  —  ffi(rm2— r2m)(ni—  n)  y™o«i n2  f0fif2 

(193)’ 

where  we  have  introduced 


K w!bu Pi  =  ml,  m!b  £  Z,  (186) 

where 

wo  =  (l)  -P1P2/K)  , 

iu[  =  (0,1,0),  «,'  =  (0,0,1).  (187) 

Inverting  this  relation,  the  eigenvalues  are  given  by 

K/31  =  M0  -|(QxM), 

Kd]  =  M2,  K^  =  -Ml5 
K/302  =  Q0-|(Mxg), 

K/32  =  Q2,  K/32  =  -Qi,  (188) 

where  M„  and  are  integers. 


B.  The  surface  partition  functions 

With  the  twist  (166),  the  two  copies  of  the  surface 
theories  are  coupled  together.  The  partition  functions 
are  given  by 


ynon,n2  ’7r0r1r2 
/Jr\r2  n nm2 


(189) 


ho  =  n0  —  r2t\  +  r2Si  mod  K, 
h 0  =  r0  +  n2t\  —  n2s\  mod  K.  (194) 

For  the  above  equation,  if  we  write  down  the  phase  in 
terms  of  ho  and  fo,  it  will  be  independent  of  t,  and  s*. 
In  other  words,  for  two  different  Xn°nin2 Xr°rir2 ,  if  they 
have  the  same  h 0,  ho,  n±,  r  1,  n2  and  r2,  the  phases  they 
acquire  under  the  U2  transformation  are  the  same.  This 
motivates  us  to  combine  these  partition  functions  and 
define,  for  fixed  ho,  ho  £  Zk, 


v. 

/\-T 


E  ^ 

*l,2>Sl,2^^K 


^■n0n1n2  j^ronr2 


(195) 


where  the  sum  is  taken  over  all  quartets  (ti,2,Si)2)  giv¬ 
ing  rise  to  given  ho, ho-  Observe  that  Xn°nin2  is  labeled 
by  two  Zk  x  Zi<-valued  quantum  numbers,  and  one  Zk- 
valued  quantum  number.  On  the  other  hand,  Xf^V"2 
depends  on  six  Zj<  x  Zx-valued  indices.  There  are  K6 
sectors.  From  Eq.  (193),  it  is  straightforward  to  read  off 
the  transformation  of  x  under  the  U2  transformation: 

(^X)£“ 

_  g—  {noni+rori)  —  ^£  (rin2  —  r2ni)(ni  —  ri)^nomn2  (  ^90) 


where,  as  before,  we  decompose  the  quantum  numbers  as 

M/l  =  KNti+nfl,  =  0, 1, . . . ,  K  —  1,  1VM  £  Z, 

Qn  =  +  rM,  rM  =  0, 1, . . .  ,K  -  1,  £  Z, 

(190) 

and  noted,  following  the  discussion  in  Sec.  IV  B  2,  the 
partition  functions  depend  only  on  the  fractional  parts 
of  M^/K  and  Q^/K.  Following  Sec.  IV  B  2  further,  we 
can  write  the  partition  functions  as 


As  for  the  U[  transformation,  the  product 
xn0n1n2Xf0f1f2  transforms  under  U[  as 

( 'JJ^X^n0nin2  (JJ^xy0^1^2 


where  the  phase  9  is  given  by 
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fl  xnon'in2  xror'ir2  (197) 


urA  iGZk 


yn0nin2  yr*orir2 
Zjr\V2  Zjn\ri2 
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(191) 


where 


n0  :=  no  +  s  x  n  +  r  x  t  mod  K, 
r0  :=  ro  +  t  x  r  +  n  x  s  mod  K, 


no  =  Kn0  +  (r  x  n) 
hi  =  Ktr  +  ni 

h2  =  K  t2  +  n2 


f0  =  Kr0  +  (n  x  r) 
fi  =  Ksi  +  ri 

f2  =  Ks2  +  r2 


(192) 


flfi'slf'fi  =  2?mo<  ,  27min'o  ,  27rAV'i  ,  2irrir'0 

K2  K2  K2  K2 

27T 

=  -^*(n0n'i  +  h'0ni  +  h0r[  +  f^ri) 

O'TT 

+  j^[(r  x  n)(n[  -  r[) 

+  ( r[n2  -  r2n'i){m  -  n)] .  (198) 

To  derive  this  result,  we  note  that  n2  and  r2  are  invariant 
under  the  U[  transformation.  As  in  our  previous  discus¬ 
sion  on  the  U2  transformation,  it  is  crucial  to  observe  that 
the  phase  9n onirori  is  independent  of  ti>2,si,2.  We  are 
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(n'1,n;n2) 

(n[,r1;r2) 

(r'i,ri;n2) 

(ni,m;r2) 

27rno  n\  r-\ 

27T ro.n\  ri 

A-nnor'-t  ri 

Anron'-,  m 

K2 

K2 

K2 

K2 

(ni,n;n2) 

(ni,n;r2) 

(r(,n;n2) 

(ni,m;r2) 

27rci2n2riirl 

27rqir2?i^T’i 

47rqi  n2‘r'i  r\ 

47rq27"2ni  nl 

K2 

K2 

K2 

K2 

TABLE  I.  The  braiding  statistical  phases  (the  second  line) 
for  the  braiding  processes  between  loop  a  and  loop  b  with 
base  loop  c  linking  both  of  them  (denoted  by  (a,  b;  c)  in  the 
first  line).  Here,  a,  b  and  c  are  the  quantum  numbers  for  loop 
excitations. 


thus  led  to  consider  the  partition  functions  x””"1,."2  de- 
fined  in  Eq.  (195),  which  transform,  under  the  U[  trans¬ 
formation,  as 


(tv 

V ^  lX-J  roriV2 


|r|2K2 


rV"l,2>I'o>rl,2 


v  y  7  72  A  ,  A 

Af'r  \r'  un1,n2ur1,r2 


(199) 


Summarizing,  the  modular  S  and  T  matrices  are  given 
by 


Sn 


—  A  ,  A  .  p-2Ki(fto"2-fio™'i+f;or2-i:or() 

un1,n2urllr2c 
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-  {n0n1+f0r1)+j^(r1n2-r2n1)(n1-r1) 


(200) 


where  n ^  =  (no,  111,712).  Observe  that  a  quick  way  to 
obtain  this  three  loop  braiding  phase  is  to  replace  no  — > 
no  +  (r  x  n)/K,  ro  — t  ro  +  (n  x  r)/K,  in  the  S  and  T 
matrices  for  the  surface  of  the  BF  theory  (Eq.  (113)). 

The  first  exponential  in  the  S  matrix, 
e-2TTx(n'0n2  )/k^  anc[  ^he  first  term  in  Eq.  (198) 
represents  the  particle-loop  braiding  phase,  which 
exists  also  in  the  ordinary  BF  theory.  On  the 
other  hand,  the  second  exponential  in  the  S  ma¬ 
trix,  e-^H(n1+r1)(n2r'1+n[r2)  —  -]/K2^  &nd  the  seCOnd 
term  in  Eq.  (198)  describes  a  topological  invariant 
which  can  be  considered  as  the  higher  dimensional 
generalization  of  the  linking  number  of  closed  lines  (in 
three  dimensions),  and  is  also  related  with  the  three- loop 
braiding  process.34-3'  More  precisely,  from  the  second 
term  in  Eq.  (198),  one  can  extract  three-loop  braid¬ 
ing  statistical  phases.  For  example,  the  phase  factor 
e2-Kir1n2n1/K  in(duded  jn  Eq.  (198)  can  be  interpreted 
as  the  three-loop  braiding  statistical  phases  associated 
to  two  loops  running  in  the  x-direction  with  quantum 
numbers  r\  and  n\  in  the  presence  of  a  base  loop 
running  in  the  y-direction  with  quantum  number  712 
(Table  I).  The  three-loop  braiding  statistics  encoded 
in  the  5-matrix  can  be  further  understood  through 
dimensional  reduction  discussed  below. 

As  for  the  T  matrix,  the  first  phase  factor 
e27rz(n0ni+r0ri)/K  pr0p0sed  to  be  the  topological  spin 
for  the  composite  particle-loop  excitations  in  the  BF 


TABLE  II.  Same  as  Table  I,  but  for  generic  values  of  the 
parameter  qij2  =  0, 1,  2, . . . ,  K  —  1. 


theory.  On  the  other  hand,  the  second  phase  factor 
e2m(r1n2-r2n1)(n1-r1)/K2  can  de  considered  as  the  topo¬ 
logical  spin  for  the  loop  excitations  with  a  base  loop 
threading  through  it.  For  instance,  e2,r!ri"2"i/K"  rep¬ 
resents  the  topological  spin  for  the  loop  excitation  with 
quantum  number  (ri,  ni)  threaded  by  the  loop  excitation 
carrying  quantum  number  77.2  ■ 

These  results  extracted  from  the  boundary  S  and  T 
matrices,  (200),  are  consistent  with  the  previous  bulk 
calculations  in  the  literature. 34-37  In  particular,  in  Ref. 
34  the  S  and  T  matrices  in  the  bulk  are  calculated  in 
the  basis  that  is  constructed  from  the  so-called  minimum 
entropy  states  (MESs)  on  the  bulk  spatial  three  torus.  In 
Ref.  36,  the  bulk  S  and  T  matrices  were  constructed  for 
.ZAq  x  Z;y2  x  Zn3  gauge  theories. 

Several  comments  are  in  order: 

(i)  The  entropy  St  computed  from  these  characters 

and  the  modular  S  matrix  shows,  in  the  limit  Ri/Rq  — > 
00  and  R1/R2  00,  the  asymptotic  behavior  St  = 

— 2  In  K  +  •  •  • ,  where  •  •  •  is  the  term  proportional  to  the 
area  of  the  surface.  I.e.,  the  constant  piece  in  the  (entan¬ 
glement)  entropy  is  the  same  as  the  two  decoupled  copies 
of  the  BF  theories. 

(ii)  For  (3  +  l)d  topological  phases  (gauge  theories) 
with  Zk  x  Zk  gauge  symmetry,  we  expect  there  are 
(at  least)  K2  different  topological  phases  that  are  differ¬ 
entiated  by  their  three-loop  braiding  statistical  phases. 
This  is  expected  from  the  group  cohomology  classification 
(construction)  of  SPT  phases;  from  U4[ Zk  x  Zk,  C/(l)]  = 
Zk  x  Zk,  we  expect  there  are  at  least  K2  different  SPT 
phases  in  (3+l)d  protected  by  unitary  on-site  symmetry 
G  =  Zk  x  Zk-  Once  the  global  symmetry  is  gauged, 
these  different  SPT  phases  give  rise  to  K2  different  topo¬ 
logically  ordered  phases  which  are  differentiated  by  the 
three-loop  braiding  phases.58,59  The  model  we  studied  in 
this  section,  the  two  copies  of  coupled  BF  surface  theo¬ 
ries,  corresponds  to  the  surface  theory  of  one  of  the  K2 
topological  phases.  The  surface  theories  of  all  the  other 
topological  phases  can  be  obtained  by  tuning  the  coeffi¬ 
cient  in  the  coupling  terms.  In  our  model,  the  coefficient 
Pi  and  P2  in  front  of  the  cubic  terms  are  chosen  to  be 
K.  In  general,  they  can  take  value  qiK  and  q2K  with 
q12  =  0, 1, ...  ,K  —  l30  which  will  lead  to  K2  different 
topological  phases  with  different  S  and  T  matrices.  The 
three-loop  braiding  phases  will  be  slightly  modified  and 
are  shown  in  Table  II,  which  are  consistent  with  Ref.  35. 

Observe  also  that  for  G  =  Zk,  IF4[Zk,  17(1)]  =  0,  i.e., 
there  is  no  non-trivial  SPT  phase  protected  by  G  =  Zk 
symmetry.  Hence,  there  is  essentially  only  one  topologi¬ 
cally  ordered  phase  with  Zk  gauge  group,  whose  surface 
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is  described  by  the  one-component  surface  theory  studied 
in  Sec.  III.  On  the  other  hand,  the  two-component  surface 
theory  studied  in  this  section  allows  richer  possibilities. 

(iii)  An  insight  on  the  three-loop  braiding  statistics 
phase  can  be  obtained  from  dimensional  reduction.  For 
the  trivial  two-component  BF  theory,  there  is  only  a  non¬ 
trivial  particle  and  loop  braiding  phase  described  in  Eq. 
(113).  This  model,  after  dimensional  reduction,  reduces 
to  the  D(Z k  x  Zk)  quantum  double  model  with  the  K- 
matrix  given  by  Kcr^  ®  Ker^, . 

For  the  topological  phase  with  non-trivial  three-loop 
braiding  statistics  phase,  the  dimensional  reduction  is 
more  interesting.  Here,  we  consider  the  simplest  non¬ 
trivial  example  with  K  =  2.  We  perform  dimensional  re¬ 
duction  on  U[  and  T  defined  in  Eq.  (199)  and  Eq.  (200). 
When  we  do  so,  we  need  to  fix  the  quantum  numbers 
n2  and  r2.  For  example,  if  we  take  ri2  =  0  and  r2  =  0, 
i.e. ,  there  is  no  third  loop  connecting  the  first  and  second 
loops,  the  S  and  T  matrices  after  dimensional  reduction 
are  the  same  as  those  for  the  two  copies  of  the  toric  code 
model. 

On  the  other  hand,  if  we  take  ri2  =  1  and  r2  =  1,  the 
dimensional  reduction  results  in  the  S  and  T  matrices 
given  by 


^rii, ri 
/  rii  ,ri 


—  ^f,Tri(non'1+n'0ni+r0r'1+r'0r1)+iri(r1-n1)(nl1  —  r'1) 

~  4 

=  5niyAiyeni(noni+rori)-^(m~ri)2-  (201) 


This  indicates  that  the  (2+l)d  topological  order  de¬ 
scribed  by  the  K-matrix 


/  2  2  -2  0\ 
2  0  0  0 
-2  0  2  2 
\  0  0  2  0/ 


(202) 


By  an  SL( 2,Z)  similarity  transformation,  this  K-matrix 
is  equivalent  to  K  =  2 <rz  ©  2<t2,  which  represents  two 
copies  of  the  double  semion  model. 

Similarly,  if  we  choose  (n2,r2)  =  (1,0)  and  (0,1),  the 
corresponding  (2  +  l)d  topological  order  is  described  by 
the  K-matrix 


and 


/  0  2  -1  o\ 
2  0  0  0 
-10  2  2 
\  0  0  2  0/ 


/  2  2  -1  0\ 
2  0  0  0 
-10  0  2 
\  0  0  2  0/ 


(203) 


(204) 


respectively.  To  summarize,  after  dimensional  reduction, 
the  original  (3  +  l)d  topological  order  with  non-trivial 


three-loop  braiding  statistics  “splits”  into  four  differ¬ 
ent  (2  +  l)d  topological  order,  which  are  controlled  by 
the  quantum  numbers  n2  and  r2.  This  result  seems  to 
be  related  with  the  group  cohomology  classification  of 
symmetry-protected  topological  (SPT)  phases  in  (2  +  l)d 
with  G  =  Z2  x  Z2  symmetry,  i.e.,  H3[ Z2  x  Z2,f/(1)]  = 
Z2  x  Z2. 


VI.  DISCUSSION 


Let  us  summarize  our  main  results. 

-  In  the  (3+l)d  BF  theory,  we  have  demonstrated, 
through  explicit  calculations  in  the  boundary  field  theo¬ 
ries  and  by  comparisons  with  known  bulk  results,  there 
is  a  bulk-boundary  correspondence  in  (3+l)d  topological 
phases.  In  particular  the  modular  S  and  T  matrices  are 
calculated  from  the  gapless  boundary  field  theory  and 
shown  to  match  with  the  bulk  results. 

-  The  surface  theory  of  the  (3+l)d  BF  theory  with 
the  theta  term  is  introduced  and  solved.  The  action  of 
the  modular  S  and  T  transformations  on  the  partition 
functions  is  calculated.  It  is  shown  that  the  partition 
functions  do  not  form  the  complete  basis  under  the  mod¬ 
ular  S  and  T  transformations. 

-  Finally,  we  propose  a  (3+1  )d  bulk  field  theory 
with  cubic  coupling  that  may  realize  three-loop  braid¬ 
ing  statistics.  We  discuss  the  twist  that  the  cubic  term 
of  the  field  theory  adds  to  the  zero  modes.  By  consider¬ 
ing  the  alternative  form  of  the  bulk  and  boundary  field 
theories,  in  which  the  quantization  rule  of  the  zero  modes 
is  twisted,  we  computed  the  surface  partition  functions, 
and  the  S  and  T  matrices  are  constructed. 

These  results  extend  the  well-established  bulk¬ 
boundary  correspondence  in  (2+l)d  topological  phases 
and  their  (l+l)d  edge  theories.  Our  approach  from  the 
surface  field  theories  provide  an  alternative  point  of  view 
to  (3+l)d  topological  phases,  and  to  recently  discussed, 
novel  braiding  properties,  such  as  three-loop  braiding 
statistics. 

There  are,  however,  still  some  aspects  in  the  (2+l)d- 
(l+l)d  correspondence,  which  we  do  not  know  if  have  an 
analogue  in  the  (3+l)d-(2+l)d  correspondence.  For  ex¬ 
ample,  in  the  case  of  the  bulk-boundary  correspondence 
connecting  (2+l)d  topological  phases  and  (l+l)d  edge 
theories,  that  the  edge  theories  are  invariant  under  an 
infinite-dimensional  algebra  seems  to  play  a  significant 
role:  the  Virasoro  algebra  or  an  extended  chiral  algebra 
of  (l+l)d  CFTs  faithfully  mirrors  bulk  topological  prop¬ 
erties  of  (2+l)d  bulk  phases.  On  the  practical  side,  that 
edge  theories  enjoy  an  infinite-dimensional  symmetry  al¬ 
gebra  provides  many  non-trivial  solvable  examples.  For 
our  example  of  (2+l)d  surface  theories  of  (3+l)d  topo¬ 
logical  phases,  on  the  other  hand,  we  did  not  make  use 
of  such  infinite-dimensional  symmetry.  In  fact,  the  sur¬ 
face  theories  studied  in  this  paper  are  not  conformal  field 
theories.  For  example,  the  two-point  correlation  function 
of  the  boson  field  (c p(t ,  r)<j>(t',  r'))  in  the  free  boson  the- 
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ory  in  (2+l)d  decays  algebraically.  This  should  be  con¬ 
trasted  with  the  logarithmic  decay  of  the  corresponding 
correlator  in  the  (l+l)d  compactified  boson  theory.  As  a 
consequence,  the  correlation  functions  of  the  bosonic  ex¬ 
ponents  exp[zm^(t,  r)]  (to  €  Z)  do  not  decay  algebraically 
in  the  (3+l)d  free  boson  theory.  Whether  or  not  there 
exists  a  unified  field  theory  framework  in  (2+l)d  field 
theories  that  strongly  resonates  with  topological  proper¬ 
ties  of  (3+l)d  bulk  topological  phases  requires  further 
investigations. 
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Appendix  A:  The  surface  theory  of  the  BF  theory 
with  theta  term 


In  this  appendix,  we  go  through  canonical  analysis  of 
the  surface  theory  of  the  BF  theory  with  theta  term.  It 
is  described  by  the  Lagrangian  density 


C  = 


2n  e*j  Cj  +  ^2  e*j diQdj  V5 


dtp 


1 

2Ai 


(eijdiQ)2  -  ^—G^dupdjip.  (Al) 

ZA2 


In  this  theory,  physical  observables  are  bosonic  exponents 
exp[irmp(t,  r)],  and  Wilson  loops  expim  fc  dxlQ(t,r), 
where  m  £  Z,  and  C  is  a  closed  loop.  The  boson  Held  ip 
and  gauge  fields  Q  are  compactified  accordingly. 

The  canonical  commutators  are 


[eijdlQ(t,  r),  q mdiCm{t,  r')]  =  ( r  -  r'), 

O'jr'j 

[<P(t,r),eijdiCj(t,r')]  =  —  d(2)(r-  r').  (A2) 


The  mode  expansions  consistent  with  the  equations  of 


motion  are  given  by  (only  the  oscillator  parts  are  shown) 


r)  = 


1 


1 


ki 


C ij  >  0  — 


-V 

V2K^X1R1R2  ^  y/A(k) 
it(k)e*A+(k)t  +  a(— k)e*A_(k)t 
Ai  v  -  1 


E 



8tt2R1R2  ^  yj A(k) 

lA+(k)at(k)eiA+(k)t  +zA_(k)a(— k)eiA-(k)t’  . 

(A3) 


where 


^2(k)  =  S(k)  =  ^eydiQkj, 


K2AiA2  ”  l 

A(k)2  =  S(k)2  +  A(k)2 

l  / 


K2Ax 


A2 


K2Ai 

A±(k)  =  -H(k)+A(k) 
where 


(2?r)E^  +  )  kikj 


K2Ax ' 


jjij  _  (  Q2/R2  —Q1Q2/R1R2 

Q1Q2/R1R2  Ql/R\  . 
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(A5) 


We  now  attempt  to  derive  the  zero  mode  quantization 
rule  solely  from  the  surface  theory.  The  canonical  com- 
muators  among  the  zero  modes  are 


-1 


[<V  PA  =  iMHu 


(A6) 


where 


f  K 

0 

0  1 
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0  1 

M  = 

-PQ2 
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K 

0 ) 

(A7) 

The  commutators  take  a  canonical  form  in  the  rotated 
basis 


j3\  =  P„Mux  =>  [a^,Px  (A8) 

Thus,  from  the  compactification  condition  on  aM,  the 
eigenvalues  of  /?  are  integers: 

Px  =  mx  S  Z.  (A9) 

This  means  the  eigenvalues  of  /3  are 

Pi*  =  P\M =  mxMXll , 

=  (l?  +  i&(miQl  +  TO2Q2),?,^i)  •  (A10) 

Renaming  integers  by  too  — >  no,  mi  — >  —  n2,  m2  —>  n\, 
(Po  ,A-«  =  g-|(Qxn),|,|).  (All) 

This  is  consistent  with  the  quantization  rules  derived 
from  the  equations  of  motion  in  the  bulk  field  theory. 
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Appendix  B:  5-function  forms 

In  this  appendix,  we  summarize  the  properties  of  6- 
function  forms.  For  an  ?r-dimensional  submanifold  A f  of 
a  D-dimensional  maniofold  A4 ,  we  define  a  ( D  —  n)-form 
6D-n(AT)  by 

[  An  =  f  So-niN)  A  An ,  VA„,  (Bl) 
J  N  Jm 

where  An  is  an  arbitrary  n-form  on  A4.  If  we  flip  the 
orientation  of  J\f , 

<^D-n(—  Af)  =  —SD-n{Af).  (B2) 

More  generally,  for  oriented  submanifolds  A/i, 

KY,c^)  =  Ec^(A/'i)  (B3) 

%  ^ 

where  Cj  is  a  coefficient. 

The  exterior  derivative  acts  on  the  delta  function  form 
as 

SD-n+1(dAf)  =  (-l)D-n+1dSD-n(AT).  (B4) 


Let  A/"i  and  A/2  be  a  submanifold  of  Ad  with  dimensions 
Tii  and  ri2,  respectively.  Define  d  as 

d  =  n\  +  ri2  —  D.  (B5) 

When  d  >  0,  A/i  and  A/2  can  have  a  d-dimensional  inter¬ 
section  within  M.  By  properly  defining  an  orientation, 
we  define  the  intersection  of  J\fi  and  A/2,  I  =  A/i#A/2- 
The  orientation  of  I  is  defined  to  be  consistent  with 

8d-<i(I )  =  £d-tu(A/i)  A  />D-ni(A/2).  (B6) 

If  the  two  submanifolds  have  complementary  dimensions, 
ni  +  ri2  =  D,  they  intersect  at  points.  Then,  the  inter¬ 
section  number 

/ (A/i,  A/2)  =  j  <5„(M)  A  8D-n(N2)  G  Z  (B7) 

counts  the  number  of  intersection  points. 

The  linking  number  of  two  submanifolds,  £  and  A/",  is 
defined  when  dim£  +  dim  A/”  +  1  =  D.  By  considering 
an  auxiliary  manifold  satisfying  <9A/"i  =  C,  the  linking 
number  is  given  by 

Lk(£,Af)  =  /(A/ijAA).  (B8) 
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